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-  i- 
Abs tract 
Hill's  equation  with  a  periodic  coefficient  of  period  n  is  investigated. 
The  problem  of  the  coexistence  of  two  linearly  independent  solutions  with  period 
n  or  2n  is  being  considered,  and  all  known  cases  in  which  such  a  coexistence 
occurs  are  shown  to  be  special  cases  of  a  four-parametric  equation  which  can  be 
transformed  into  an  equation  of  Hill's  type.  Also,  the  question  of  the 
existence  of  solutions  admitting  a  terminating  Fourier  expansion  is  being 
investigated,  and  partial  answers  are  being  provided. 
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I,  Introduction 

1.  Statement  of  the  problem 

The  general  theory  of  differential  equations  with  periodic 
coefficients  dates  back  to  Hermite  (1877)  and  Picard  (1879)  for 
doubly-perdodic  coefficients,  and  Floquet  (1883)  for  singly-period 
coefficients  (see  [26],  p. 375).  Although  a  differential  equation  with 
periodic  coefficients  does  not  of  necessity  possess  periodic  solutions, 
it  can  be  expected  that  for  some  values  of  the  parameters  of  the  equation 
such  solutions  will  be  obtained.  The  equation  of  second-order  which 
appears  to  be  most  frequently  treated  is 

(1)  y"  +  (h  +  0(x))y  «  0 

with       0(x+n)  »  0(x)  ;         0(x)dx  »  0  . 

o 

Equation  (l)  is  the  Liouville  normal  form  of  the  classical  Sturrn-Liouville 
equation  and  includes  the  well-known  equations  of  Mathieu  and  Hill« 
The  question  of  stability  of  the  solutions,  which  was  initiated  by 
Liapounov  and  Poincare,  is  of  considerable  interest.   It  turns  out  that 
there  are  infinite  sequences  of  characteristic  values  of  h 

(2)  hg  <  h^  <  h_  <  h,  <  h,  <  h^  <  h.  <  ... 

(3 )  E.  <  Ep  <  K^  <  El  <  Ej^  <  K,  <  , . , 
such  that  the  intervals 

(h)     (  -  00,  hQ),  (Ep  Ej),  (h^,  h^),  (E^,  hj^),  (h^,  h,^),  ... 


The  transformation  equations  are  given  in    [39]   and    [27] ♦ 
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are  intervals  of  instability,  and  for  values  of  h  in  these  intervals 

there  exist  solutions  of   (1)  which  are  unbounded  as     x     approaches 
infinityM.t32l^ 

The  values  of     h     given  in  (2)  are  obtained  from  the  standard 
solutions  y,  (x,h)   and  l/2(x,h)  vrtiich  are  defined  by  the  initial  conditions 

y^(0,h)  -  1  y^(0,h)  -  0 

y2(o,h)  -  0  y2(o,h)  -  0 

by  solving  the  equation 

(5)  yi(n,b)  +  7^  (n,h)  -  2  -  0  , 
while  those  of  (3)  are  obtained  from 

(6)  y^^Cn.h)  ♦  y^  («,h)  +  2-0 

and  both  sequences  are  connected  by  the  inequalities 

(7)  -  00  <  hg  <  K.  <  H"p  <  h-  <  hp  <  K  <  K,  <  h-  <  h,  <  .., 

Thus  the  solutions  with  the  fundamental  periods  n  and  2n  play  a  special 
role  in  the  stability  question  since  the  h  of  (2)  are  the  characteristic 
valiies  for  the  solutions  of  period  n,  and  the  E   of  (3)  are  the 
characteristic  values  for  the  solutions  of  period  2n.  If  the  intervals 
of  (k)  are  numbered  0,  1,  2,  3>  U,  etc.,  then  the  boundaries  of  the  even 
intervals  of  instability  correspond  to  solutions  of  period  n  while  the 
boundaries  of  the  odd  intervals  correspond  to  solutions  of  period  2n. 

The  present  investigation  arose  from  the  following  two  questions: 

(i)  Are  there  other  equations  besides  the  Lame  equation  for  which 
all  but  a  finite  nvimber  of  intervals  of  instability  vanish. 

(ii)  For  equations  in  which  0(x)  is  a  finite  polynomial  in  cos  2x 
(e.g.,  the  Mathieu  equation),  will  any  of  the  intervals  of  instability 
vanish. 
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These  questions  are  answered  by  considering  the  simultaneous 

occurrence  of  two  linearly  independent  solutions  of  period  n  or  2h.  Such 

solutions  occur  whenever  two  characteristic  values  of  h  coalesce  and 

are  always  even  and  odd  solutions  of  the  same  period.  Details  of  the 

existing  theory  and  previous  work  may  be  found  in  the  recent  review  of 

Hill's  equation  by  Magnus  and  ShenitzerL  -' , 

2.  Review  of  previous  work 

Perhaps  the  best  known  result  in  this  field  is  the  impossibility 
of  the  coexistence  of  two  Mathieu  fimctions;  that  is,  two  linearly 


independent  solutions  of  the  Mathieu  •  equation  with  period  n  or   2n  will 

not  coexist.  This  was  first  proved  by  Ince'-  -•  and  also  subsequently  by 

Mr33l 
w^^^^     and  Markovic*-  -• .  This  result  may  have  been  somewhat  surprising 

since  only  a  vear  previously  Poole'-  -■  demonstrated  the  simultaneous 

existence  of  solutions  of  period  bn  for  the  Mathieu  equation.  This  can 

be  extended  to  show  coexistence  of  solutions  with  period  mn,  m>2j  in 

Pjll 
fact,  it  is  easily  shown^  -•  that  if  Mathieu 's  equation  has  a  periodic 

solution  with  a  period  ran,  m>2,  then  all  of  its  solutions  have  this 

period. 

Undoubtedly  the  impossibility  of  coexistence  for  the  Mathieu 
equation  led  Ince  astray,  for  in  1926  InceL^  -' »  L^^-l  followed  somewhat 
later  by  Markovic^  J,  published  the  erroneous  statement  that  equation  (1.1), 
with  0(x)  an  even  function,  cannot  have  more  than  one  solution  of  period 
n     or  2n,  This  error  was  noted  by  Borg  (see  [2],  p.90)  in  his  remarkable 
dissertation  where  he  gives  the  Meissner  equation  as  a  counter-example 
to  disprove  the  statement  of  Ince  and  Markovic, 

In  the  same  thesis  Borg  proved  the  following  deep  theorems  for 
equation  (i.l): 


See  also  [26]  where  the  same  result  is  given  for  the  more  general 
Hill's  equation. 
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Theorem  2»1 

A  necessary  and  sufficient  condition  that  all  odd  interv:al3 
of  instability  shall  vanish  is  that 
0(x  +  n/2)  -  0(x) 

Theorem  2«2 

All  intervals  of  instability,  except  the  zero-th,  vanish  only 
if  0(x)  is  a  constant. 

It  is  somewhat  ironical  that  Ince*  furnished  a  counter-example  to 
his  1926  statement  in  his  second  paper  on  Periodic  Lame  Functions. 
Although  he  did  not  state  it  explicitly,  his  results  show  that,  for 
suitable  values  of  the  parameters  in  the  Lame'  equation,  it  is  possible 
for  all  but  a  finite  number  of  intervals  of  instability  to  vanish.  Since 

the  Lame"  equation  has  doubly- per iodic  coefficients,  solutions  with 

'  [6] 

imaginary  periods  can  occur  and  can,  as  was  shown  by  Erdelyi^-*  coexist. 

The  stability  question  for  the  Whittaker  equation  was  subjected  to 
an  extensive  numerical  attack  by  Klotter  and  Kotowski**.  They  give 
stability  curves  for  a  range  of  values.  However,  for  the  range  of  values 
they  examine,  only  one  intersection  of  an  odd  and  even  solution  of  the 
same  fundamental  period  (indicating  coexistence)  is  exhibited. 

Another  equation,  which  McLachlan  ^  -^   calls  the  Frequency 

..«.^uLj.-v^^..  ^^^ — ^^.., J   •  For  this  equation 

all  even  intervals  of  instability,  except  the  zero-th,  can  vanish. 
However,  none  of  the  odd  intervals  will  vanish.  These  results  are 


See  [2liJ,  in  particular,  paragraph  5,  pp.  87-89» 
See  [28j.  The  fundamental  periods  in  this  paper  are  2n  and  Un« 
Klotter  and  Kotowski,  [28],  refer  to  an  unpublished  report  of 
A.  Weigand  on  the  stability  question  for  this  equation. 
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valid  as  long  as  the  equation  remains  non-singular.   It  may  be  noted  that 
if  the  absolute  value  of  the  modulation  parameter  is  equal  to  or  greater 
than  unity,  the  equation  becomes  singular.  In  this  sense,  the  Frequency 
Modulation  equation  belongs  to  a  different  category  than  the  Mathieu 
or  Whittaker  equations. 

3.  The  classification  of  equations 

Equation  (1.1)  is  called  Hill's  equation  by  many  writers.  We  shall 
use  the  name.  Hill's  equation,  in  a  more  restricted  sense  to  designate 
the  equation 

oo      „ 

(1)  y"  +  (  £^  A  cos'^^  x)y  -  0  . 

o 

This  avoids  the  cumbersome  terminology  of  calling  (1)  the  Even  Hill's 
equation  or  the  Restricted  Hill's  equation.  If,  in  place  of  the 
infinite  series,  the  coefficient  of  y  is  a  finite  polynomial  with  m 
terms  in  cos  2x,  then  the  equation  will  be  called  the  Finite  Hill 
equation.  When  m»l  it  becomes  the  Kathieu  equation,  and  for  m«2 
it  beccmies  the  Whittaker  equation.  Another  equation  with  periodic 
coefficients  (doubly-periodic,  in  fact)  is  the  Generalized  Lame  equation 

(2)  y»  ♦  (  r"  A  sn  "  x)v  -  0  , 

o   " 

where  sn  is  the  elliptic  sine  function.  For  m"l,  equation  (2 )  reduces 
to  the  Lame  equation,  while  for  m»2,  it  is  the  Lame  wave  equation. 
An  equation  of  considerable  imoortance  in  this  study  is 

(3)  (1  +  a  cos  2x)y"  ♦  (b  sin  2x)y'  +  (c  +  d  cos  2x)y  •  0 
which  does  not  appear  to  have  been  named  previously.  We  propose  to  name 
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it  the  Ince  equation  in  honor  of  the  British  mathematician  who  contrib- 
uted so  much  to  this  field.  The  Ince  equation  is  the  most  general 
trigonometric  equation  leading  directly  without  transformation  to 
three-term  recurrence  relationships. 

It  is,  of  course,  possible  to  continue  and  list  all  the  equations 
to  be  treated  in  this  thesis.  A  more  useful  approach  is  to  display 
the  equations  in  a  sequence  of  three-dimensional  manifolds  in  a  manner 
which  indicates  their  relation  to  one  another  and  at  the  same  time  their 
position  in  the  general  scheme  of  differential  equations.  Such  a 
presentation  of  the  various  equations  may  be  found  in  the  Appendix. 

U.  Remarks  on  the  method 

The  method  invented  by  Ince''  -^  .to  show  the  impossibility  of 
coexistence  for  the  Mathieu  equation  is  simple  and  direct.  It  consists 
essentially  of  assuming  the  simultaneous  existence  of  an  even  and  an  odd 
solution  with  a  fundamental  period.  If  these  solutions  are  expressed  as 
infinite  series  of  cosines  and  sines, respectively,  the  convergence 
requirement  leads  to  the  condition  that  the  2  by  2  det/enninant  formed 
from  the  nth  and  (n  -  l)th  coefficient  of  each  seirLes,  as  n 
approaches  oo,  must  vanish.  Each  successive  determinant  is  connected 
by  the  recurrence  relationships  obtained  by  substituting  the  assumed 
solutions  into  the  equation.  From  these  relationships  it  can  be  shown 
in  the  case  of  the  Mathieu  equation  that  the  condition  for  the  vanishing 
of  the  determinant  is  that  either  the  sine  or  the  cosine  series  vanishs 
identically.  Hence,  coexistence  cannot  occur.  The  same  method  for  the 
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Lame  equation*  reveals  the  conditions  necessary  for  coexistence.  The 
success  of  the  method  depends  in  each  case  on  the  occurrence  of 
three-term  recurrence  relationships  which  lead  to  the  simple  but 
necessary  connection  between  successive  determinants. 

It  seemed  natural  and  desirable  to  exploit  this  scheme.  The 
procedure  adopted  is  to  transform  an  equation  to  sane  form  of  the  Ince 
equation  and  then  apply  the  Ince  method.  Transformation  of  the  depen- 
dent or  independent  variable,  or  both, may  be  employed,  Transfonnations 
of  the  dependent  variable  must,  of  course,  be  periodic. 

A  possible  extension  of  the  Ince  method  to  the  case  in  which 
four-term  recurrence  relationships  occur  is  introduced.  The  idea  is  to 
use  3  by  3  determinants  in  place  of  2  by  2  determinants.  The  full 
scope  and  validity  of  this  extension  has  not  yet  been  determined. 

S,     Scope  of  method 

In  this  paper  the  two  questions  of  Section  1  are  answered 
affirmatively.  The  second  question  is  answered  by  the  discussion  of  the 
Whittaker  equation  which  is  next  after  the  Kathieu  equation  in  the 
sequence  of  Finite  Hill  equations.  The  Whittaker  equation  would  normally 
lead  to  five-term  recurrence  relationships'-  -',  and  hence  it  is  trans- 
formed to  an  equation  yielding  three-term  recurrence  relationships.  The 
necessary  and  sufficient  conditions  for  coexistence  of  solutions  of  the 
Whittaker  equation  are  developed  from  this  transformed  equation  which  is 
a  reduced  form  of  the  Ince  equation. 

The  Ince  equation,  which,  as  mentioned  before,  is  the  most 
general  trigonometric  equation  leading  directly  to  three-term  recurrence 

See  [2^ .  The  Lame  equation  must  first  be  transformed  into  a  special 
case  of  the  Ince  equation. 


relationshipSj provides  the  answer  to  the  first  question.  All  of  the 
equations  treated  previously,  including  the  Whittaker  equation,  contain 
two  or  three  arbitrary  p)araineters.  The  complete  Ince  equation  has 
four  arbitraiy  parameters  and  thus  represents  an  important  extension. 
There  is  also  the  possibility  of  treating  certain  equations  with  more 
than  four  parameters  if  they  can  be  transformed  to  the  Ince  equation  and 
the  excess  parameters  imbedded  in  the  transformations.  Equations  with 
Jacobian  elliptic  functions  as  coefficients  such  as  the  Lame  and  the 
Hermite  Elliptic  equation*  are  handled  by  transforming  them  into  Ince 
equations. 

Some  results  on  the  coexistence  question  are  obtained  for  the 
general  Finite  Hill  equation.  These  give  an  indication  of  the  nature  of 
the  solutions  connected  with  coexistence.  The  Finite  Hill  equation 
finds  an  important  application  in  the  study  of  non-linear  oscillations, 
since  the  question  of  stability  of  periodic  solutions  for  non-linear 
differential  equations  reduces  to  a  question  of  the  stability  of  a  Hill's 
equation  (See  [li2],  p.  192). 

Certain  singular  equations  are  discussed  briefly  and  are 
interesting  because  they  are  the  equations  which  arise  from  the  separation 
of  the  wave  equation  in  orthogonal  coordinate  systems.  An  inherent 
difficulty  is  the  problem  of  the  existence  of  even  a  single  periodic 
solution,  since,  as  is  well  known,  a  differential  equation  with  periodic 
coefficients  does  not  of  necessity  have  periodic  solutions.  The  point  of 


See  [12] .  The  name  Hermite  Elliptic  equation  was  chosen  to  avoid 
confusion  with  the  better  known  differential  equation  for  the 
classical  Hermite  polynomials. 
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view  taken  here  is  to  examine  the  possibility  of  coexistence,  asstaning 
the  existence  of  oeriodic  solutions.     Singular  equations  can  lead  to 
two-term  and  four-term  recurrence  relationships  as  well  as  three-term. 
Non-singular  equations  must  have  recurrence  relationships  with  an  odd 
number  of  terms. 
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II»  The  Whittaker  Equation 

6.  Introductory  remarks  and  origin  of  the  equation 

This  equation  was  introduced  by  E.T.  V/hittaker  in  his  1915 
paper  ^-' * '--'  on  differential  equations  whose  solutions  satisfy 
homogeneous  integral  equations.  Actually  the  equation  appeared  earlier, 
in  1902,  in  a  memoir  by  A,  Liapounof  f  L  -J,  as  an  illustrative  example 
for  the  general  theory  developed  theire.  As  the  standard  form  of  the 
Whittaker  equation,  we  shall  adopt 

2 
(1)     y"  +  (h  +  hmq  cos  2x  +  2q  cos  Ux)y  *  0 

which  (with  q  replaced  by  -p/U)  is  the  equation  given  by  Whittaker. 
There  is  some  loss  of  generality  in  this  form  for  equation  (1),  since 
the  coefficient  of  cos  Ux  is  constrained  to  be  positive.  The 
question  of  a  negative  coefficient  will  be  examined  in  Section  10. 

If  q  goes  to  zero  as  ra  approaches  infinity,  and  at  the  same 
time  raq  remains  finite,  then  (1)  becomes  the  Mathieu  equation. 
Whittaker  remarked  that  the  functions  which  satisfy  equation  (1)  are  re- 
lated to  the  Mathieu  functions  in  the  same  manner  that  the  Associated 
Legendre  functions  are  related  to  the  Bessel  functions.  Ince  pointed 
out  that  the  Whittaker  equation  has  the  same  relation  to  the  confluent 
hypergeometric  equation  as  the  Mathieu  equation  has  to  the  Bessel 

ri9i 

equation'-  -*  •  Additional  properties  of  the  Whittaker  equation  were 
discussed  by  InceL^J  in  his  paper  on  the  real  zeros  of  its  solution. 

The  Mathieu  equation,  as  is  well  known,  arose  from  the  discussion 
of  the  vibrations  in  an  elliptic  membrane  (see  [353iand  Qi^,  ch.  19). 
Therefore  it  is  not  surprising  that  the  Whittaker  equation  occurs  in  the 
physical  problem  of  the  vibrations  of  an  elliptic  membrane  with  a  square 
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r2ii 

law  decrease  in  its  density  radially  from  the  center^  -•,  Humbert 
obtained  the  Whittaker  equation  from  a  separation  of  the  Laplace  equation 
(in  four  dimensions)  by  introducing  hypercylinders  formed  by 
three-dimensional  confocal  paraboloids.  The  Mathieu  equation  arises  in  an 
analogous  manner  in  three  dimensional  space  from  the  confocal  paraboloids. 
The  Whittaker  equation  could  also  be  obtained  from  a  separation  of  the 
three-dimensional  wave  equation  in  a  region  bounded  by  confocal 
paraboloids,* 

A  more  recent  application  of  the  Whittaker  equation  is  found  in  the 
paper  by  Chang  on  the  stability  of  periodic-beam  focusing L  J,  where  the 
stability  chart  developed  by  Klotter  and  Kotowski**  is  employed, 

?♦  Periodic  solutions 

If,  under  the  assumption  that  periodic  solutions  exist,  a 
trigonometric  series  is  substituted  into  the  Whittaker  equation 

(1)  y"  +  (h  +  Imq  cos  2x  +  2q  cos  Ux)y  ■  0  , 

the  recurrence  relationships  which  result  would  have  five  terms.  In 
order  to  obtain  a  sequence  of  three-term  recurrence  relationships  it 
is  necessary  to  consider  solutions  of  the  form 

(2)  y  =  u  exp(q  cos  2x)  . 

The  periodic  solutions  of  interest  are  those  with  period  n  or  2n, 

r32l 
These  can  be  shown"-  -■  to  be  series  of  one  of  the  following  four  types: 


„  ■■"-   I     I   IB  IB  ■■  ■  ■  I    II     .1 

See  |l6].  The  remark  about  the  wave  equation  is  credited  to  E.T. 
Copson.  , 

See  [28],  p.  155.  Chang  maKes  the  change  in  scale  required  to  change 
the  last  period  from  2n  to  n. 
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00 

(Ce)      ^  ■  L-  ^2  ^°^  ^"^ 
o 


00 

(Se)      ^  ■  EI  ^2n  ^^  ^"^ 


00 

(Co)      u  -  21  A2n+1  *^°^  (2n+l)x 
o 

00 

(So)      ^  ■  EI  ^2n+l  ^^"  (2n+l)x 
o 


since  the  recurrence  relations  for  the  coefficients  in  a  general  Fourier 
expansion  connect  only  the  coefficients  appearing  in  (Ce)  etc. 

Without  loss  of  generality  it  may  be  assumed  that  for  the  even 
solutions  Ce  and  Co 

(la)      u(0)  -  1    and    u'(0)  =  0 
while  for  the  odd  solutions  Se  and  So 

(lb)      u(0)  -  0    and    u'(0)  -  1  . 

The  values  of  the  parameters  for  which  these  solutions  exist  provide 
the  boundary  curves  in  the  h,q-plane  which  separate  the  intervals  of 
stability  and  instability.  Whenever  two  solutions  of  period  n  or  2n 
exist  simultaneously  for  the  same  values  of  the  parameters,  the 
corresponding  interval  of  instability  vanishes. 

The  differential  equation  obtained  by  employing  the  transformation 
(2)  is 

(3)       11"+  (a  sin  2x)u»  +  (b  +  c  cos  2x)u  »  0 

with 
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(U)     a  »  -Uq;   b»h+2q;   c-  U(ra-l)q  . 

The  choice  of  +q  in  (2)  is  arbitrary  and  no  restriction  is  implied 

for  a  which  may  be  positive  or  negative.  The  positive  sign  preceding 

2 
the  term  2q  cos  hx    is,  however,  essential. 

Consider  the  periodic  solutions  with  period  n.  By  direct  substi- 
tution of  the  series  (Ce)  in  the  differential  equation,  the  recurrence 
relations  are  obtained: 

(c-2a)A  +  2bA^  -  0 

,^^  (c-Ua)A^  +  2(b-U)A2  +  2cA^  »  0 

(c-6a)A^  +  2(b-16)A,  +  (c+2a)A2  =•  0 

and  the  general  term  for  n  >  2,  is  : 

(6)  [c  -(2n+2)a]A2^^2  +  2[b  '{2n)^}A2n   *  &  *^2n-2)a] Ag^.g  "  °  * 

If  a  periodic  solution  Ce(x)  exist?,  its  Fourier  series  will  be 
convergent  for  all  complex  values  of  x  since  Equation  (1)  has  no 
singularities.  From  this,  it  can  be  derived  that 

n  

(7)  lim    /(TT   -  0 
n  -■>  oo 

and,  in  particular,  that  A  ->  0  as  n  ->  oo. 

To  decide  for  which  values  of  the  parameters  a  periodic  solution 

exists  is  a  difficult  problem.   It  can  be  shown  (see  [32]  for  a  detailed 

proof)  that  for  fixed  values  of  m  and  q  (or,  alternatively,  of  a 

and  c)  there  exists  a  denumerable  infinitude  of  values  of  h  (or,  alternatively, 

of  b)  for  which  (1)  (or  (3))ha3  a  periodic  solution  of  type  Ce  or 

of  any  one  of  the  types  Co,  Se,  and  So.  The  corresponding  values  of 

h  (or  b)  are  the  roots  of  a  transcendental  equation  which  is  called  the 
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characteristic  equation.  Its  roots  are  all  real  and  bounded  from  below, 
and  their  only  point  of  accumalation  is  oo. 

The  characteristic  equation  for  b  may  be  written  in  a  form 
involving  a  continued  fraction* 

(8)    h  -   (l/2)(2a-e)(c)/(2)^  ^  (Ua.c)(2a^c)/(lb)^  ^ 
1  -  b/2^  1  -  b/U^ 


(6a-c)(ha-^c)/(U8)^  ^   ^  (2na-c)  ( [2n-2la-Kc)/[2(2n)(2n-2)3^  . 
1  -  h/e'^  1  -  h/{2nr 

For  the  series  (Se),  the  recurrence  relations  are 


(c-ha)B^^  +  2(b-U)32  «  0 

(9)  (c-6a)B^  +  2(b-16)B,  +  (c+2a)B^  -  0 

o  U         <i 

(c-8a)Bj,  +  2(b-36)B^  +  (c+Ua)B^  '  0 
and  the  general  term  for  n  >  2  is 

(10)  [c-(2n+2)a]B2^^2  *  ^  [b-(2n)^jB2„  +  [c+(2n-2)a]B2„,2  "  °- 

The  coefficients  of  the  B's  in  equation  (10)  arise  from  those  of 
the  A's  in  equation  (6)  by  omitting  the  first  one.  The  characteristic 
equation  for  b  is 

(11)  ^  -  2^  ,  2^(Ua-c)(2a>c)/(16)^  ^  (6a-c) (Ua^c)/(U8)^  ^   ^ 

1  -  h/>/  1  -  b/6^ 

*  (2na..c)[(2n-2)a-»-c]/[2(2n)(2n-2)]^  . 
1  -  b/(2n)^ 


The  continued  fraction  converges  for  finite  values  of  a,  b,  and  c. 


See  [373,  paragraph  56,  Theorem  Ul. 
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If,  in  (8),  the  first  complete  denominator  is  equated  to  zero, 
the  resulting  equation  for  b  is  exactly  equation  (11). 

The  continued  fractions  given  by  (8)  and  (11)  vdll  terminate 
for 

(12a)        c  -  2na  n  >  1 

(12b)        c  -  -  2na  n  >  0 

unless  the  corresponding  complete  denominator  vanishes  simultaneously. 
It  should  be  nested  that  equation  (8)  also  may  terminate  for  n  ■  0 
(i.e.jC  =  0),  and  n  "  1  (i.e,,c  =  2a).  When  the  continued  fraction 
terminates  then  the  equations  become  algebraic  and  there  are  at  most  n 
or  n+1  values  of  b  obtained  from  (8)  and  n-1  or  n  values  of  b 
obtained  from  (11)  accordingly  as  c  is  a  positive  or  negative 
multiple  of  a.  However,  there  are  also  an  infinite  number  of  values  for 
b  obtained  by  equating  the  appropiriate  complete  denominator  to  zero. 

The  nature  of  the  periodic  solutions  with  period  n  of  the 
Whittaker  equation  may  be  inferred  as  an  almost  immediate  consequence  of 
the  recurrence  relationships  and  characteristic  equations  for  b.  Some 
of  these  inferences  are  contained  in  the  following  two  theorems: 
Theorem  7.1 

For  the  solutions  with  period  n  corresponding  to  Ce,  and  of  the 
form 

00 

(13)        y  "  [  ZI  *2n  '^°3(2n)x]  exp(q  cos  2x)  , 

o 

we  have 

(i)  if  c  ■  2na,  n  >  0, 

then  either  the  infinite  cosine  series  in  (lb)  begins  with  the  term  in 
2n  (i.e.,  the  first  n  coefficients  vanish  identically)  and  the  values 
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of  b  are  obtained  by  equating  the  complete  n-th  denominator  of  the 
continued  fraction  (8)  to  zero, 

or  the  values  of  b  are  obtained  from  the  terminated  continued 
fraction,  and  have  at  most  n  distinct  values. 

(ii)   If  c  »  -2na,      n  >  0,         then 

either  the  cosine  series  in  (lU)  terminates  with  the  (n+l)-th  term  and 
b  has  at  most  n+1  distinct  values  obtained  from  the  terminated 
continued  fraction  (9), 

or  the  values  of  b  are  obtained  by  equating  the  complete 
(n+l)-th  denominator  of  the  continued  fraction  to  zero. 

Theorem  7*2 

For  the  solutions  with  period  n,  corresponding  to  Se, 

00 

(lU)    y  -  [  ^  ^2^  sin(2n)x]  exp(q  cos  2x)  . 
(i )     If  c  ■  2na,    n  >  1,  then 

either  the  infinite  sine  series  in  (lU)  begins  with  the  term  in  2n 
(i»e«,  the  first  n-1  coefficients  vanish  identically)  and  the  values 
of  b  are  obtained  by  equating  the  complete  (n-l)-th  denominator  of 
the  continued  fraction  (11 )  to  zero, 

or  the  values  of  b  are  obtained  from  the  terminated  continued 
fraction  and  have  at  most  (n-1)  distinct  values. 

(ii)    If  c  ■  -  2na,  n  >  1,  then 

either  the  sine  series  in  (lU)  terminates  with  the  n-th  term  and 
b  has  at  most  n  distinct  values  obtained  from  the  terminated  continued 
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fraction  (12), 

or  the  values  of  b  are  obtained  by  equating  the  n-th  complete 
denominator  of  the  continued  fraction  to  zero. 

Now  consider  the  periodic  solutions  with  period  2n»  For  the 
even  solution  Co,  the  corresponding  recurrence  relationships  are 

(c-3a)A^  +  [2(b-l)  +  (c-a)]A^  -  0 

(15)  (c-5a)A^  +  2(b-9)A^  +  (c+a)A^  -  0 
(c-7a)A^  +  2(b-25)Aj  ♦  (c+3a)A^  -  0 

and  the  general  term  for  n  >  1  is 

(16)  [c-(2n+3)a]A2^^3  +  2  [b-(2nn)^  A2^^^  ♦  [c+(2n-l)a]A2^^3^  -  0  . 

As  before,       /^Pn+l  "^  ^     ^^     n  -c>  oo     for  a  periodic  solution. 
Tne  characteristic  equation  for     b     is 

(17)  b  .  1  (a-c*2)  .   (3a-c)(a^c)/6^ 

'^  1  -  b/3  +... 

+     [(2n-»-l)a-cl  [(2n-l)a-t-c]/r2(2n-l)  (2n-H)]^ 

1  -  b/(2n+l)^  +,..     . 

The  recurrence  relationships  for  the  odd  solution  are 

(c-3a)B^  +  [2(b-l)-(c-a)]B^  -  0 

(18)  (c-5a)Bc  ♦  2(b-9)B^  +  (c+a)B^  -  0 
(c-7a)B^  +  2(b-25)B^  +  (c+3a)B^  -  0 

and  the  general  term  for  n  >  1  is 

(19)  [c-(2n+3)a]B2^^3  +  2[b-(2n+l)2]B2^^^  ♦  [c+(2n-:^]B2^^^  -  0  . 
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Again     'v^I     7    ->  0     as     n  ->  cxd     for  a  periodic  solution. 
The  continued  fraction  for     b     is 

(20)       b-^(c-a.2).^3a-c)(a>c)/6^ 

'^  1  -  b/3  +  ... 

[(2n-H)a-c][(2n-l)a-t-c]/[2(2n-l)(2n->-l)3^ 
..,+  1  -  b/(2n+l)^ \        +.,.     . 

Except  for  the  leading  term,  equation  (20)  is  identical  to 
equation  (17).  These  continued  fractions  will  terminate  for 

(21a)    c  -  (2n+l)a,     n  >  1 
(21b)    c  -  -(2n+l)a,    n  >  0 

unless,  as  before,  the  corresponding  complete  denominator  vanishes 
simultaneously.  Thus  b  has  at  most  n  or  n+1  values^  or  b  is 
obtained  by  equating  the  appropriate  complete  denominator  to  zero.  If 
c«a,  that  is  n  =  0  in  (21a),  then  equations  (17)  and  (20)  are  also 
identical. 

The  solutions  Co  and  So  behave  in  a  manner  similar  to  Ce  and 
3e  and  exhibit  the  same  dependence  on  whether  c/a  is  positive  or 
negative.  However,  since  both  of  the  series  Co  and  So  go  frc»n  0  to 
00,  the  results  can  be  conveniently  summai*ized  in  a  single  theorem. 

Theorem  7.3 

For  the  solutions  with  period  2n,  corresponding  to  Co 

00 

(22)  y  ■   [2Z  *2n+l  cos(2n+l)x]   exp(q  cos  2x) 

o 

or  corresponding  to  So 

00 

(23)  y  -   m  \ri^\  sin(2n+l)xj   exp(q  cos  2x) 

o 
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(i)   If  c  =  (2n+l)a,    n  >  0,    then 

either  the  infinite  series  in  (22)  or  (23)  begin  with  the  term 
B-  ,   (i.e.,  the  first  n  coefficients  vanish  identically)  and  the 
values  of  b  are  obtained  by  equating  the  comolete  n-th  denominator 
of  the  continued  fraction  (1?)  or  (20)  to  zero, 

or  the  values  of  b  are  obtained  from  the  terminated  continued 
fraction  and  have  at  most  n  distinct  values. 

(ii)  If  c  =  -(2n+l)a,   n  >  0,    then 

either  the  trigonometric  series  in  (22)  or  (23)  terminates  with  the 
(n+l)-th  tenn  and  b  has  at  most  n+1  distinct  values  obtained  from 
the  terminated  continued  fraction  (l-'^)  or  (21) 

or  the  values  of  b  are  obtained  by  equating  the  complete  (n+l)-th 
denominator  of  the  continued  fraction  to  zero, 

8.  Necessary  and  sufficient  conditions  for  coexistence  of  periodic  solutions 

The  question  of  the  simultaneous  existence  of  fundamental 
periodic  solutions  of  the  Whittaker  equation  can  be  investigated  by  the 
method  developed  by  Ince  for  the  Mathieu  equation.  This  method  stems 
from  the  observation  that  the  general  term  of  the  recurrence  relationships 
for  Ce  is  the  same  (with  B  written  for  A)  as  for  Se,  and  that  the 
one  for  Co  is  the  same  as  for  So.  Rewrite  the  general  tern  of  the 
recurrence  relationships  for  Ce  as 

(1)  [(2n+2)a-c]A2^^2/A2n"C(2n-2)a+cjA2^.2/A2„  -  2[b-(2n)^ 
and  for  Se  as 

(2)  [(2n+2)a-c]B2„^2/^2n-^^2n-2)a*c]B2^_2/B2n"2&-(2n)2]   . 


-  20  - 


Equate  the  left-hand  sides  of  equations  (1)  and  (2)|  multiply  by  A-  B«  ; 
divide  by  2a  and  rearrange  terms  to  obtain 

(3)  C(n-l)-(c/2a)j[A2„,2B2n-^2n.2^2n]  '  "  [(n-l)*(c/2a)J  f  ^B^^.^"  B^A^^J 
or,  rewriting  in  a  convenient  notation, 


ih) 


A       A 

2n+2  *2n 


^2n-2  ^2n 


(n-1)  *   (c/2a) 
(n+1)  -  (c/^a) 


A   A 

n  *2n-2i 


\  ^2n-2| 


Equation  (U)  is  a  recurrence  relation  in  determinants.  Continuing  the 
recursion  process  and  letting 


(5)  K^  -  (n-1)  +  (c/2a)  ; 


H^  -  (n+1)  -  (c/2a) 


then,  after  recalling  that  the  first  two  terms  of  Ce  and  the  first  term  of 
Se  are  not  obtained  from  (l)  and  (2),  we  have 


(6) 


A      A 
'  2n+2   2n 


'^2n+2  ^2n 


vn-1 


n 


-  2(-l)"-^  Tr^(^r>^r^Vo  • 


Now,  firom  the  theory  of  the  Gairana  function  we  have 


N 


lim     N"'^  xfr     (1  +  -)  -  — 

N->oo  n-1  "         r  (>^) 


Therefore,  we  find  that 

/^\   -1  •   /  T  ^n  (c/a)-2 
(7)   lun  (-1)  n' 


n->oo 


A      A 
2n+2   2n 


^2n+2  ^2n 


.  1  liircZ^a) 

^  r  (c/2a)    °  2 


Now  we  know  that  the  periodic  solutions  of  our  differential  equation  will  be 
differentiable  any  number  of  times  and  that,  therefore,  for  any  positive  integer  m 

lim  n"A„  -  lim  n"B.  -  0  . 
n-c>  00       n-5>  00 

Consequently,  coexistence  of  periodic  solutions  can  be  possible  only  if 

(8)     Liirs^^  A  B,  -  0  . 
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The  case  when  a«0  (with  cj^O,   i.e.,  m->oo)  corresponds  to  the 
Mathieu  equation*  From  equation  (U)  it  is  seen  that  for  a=0  the 
coefficient  of  the  right-hand  determinant  is  1.  By  continuing  the 
recursion  process  on  the  right  and  passing  to  the  limit  on  the  left 

(9)  A  B,  m  0 

0  d 

and  hence  either  A  =  0  or  B^  ■  0.  Therefore,  either  all  the  A's  or 

o  2  ' 

all  the     B's     in  the  trigonometidc  series  vanish  identically,   indicating 

(in  accordance  with  the  well-known  result)  that  simultaneous  periodic 

solutions  can  never  occur  for  Mathieu 's  equation. 

With     8   /  0,   equation  (8)  is  satisfied  only  if     A  Bp  ■  0     except 

when     (c/2a)     is  zero  or  a  negative  integer.* 

When     A  B,,   "0     then  either  the  cosine  series  or  the  sine  series 
o  d 

vanishes  unless 

(10)  c  '  2(n+l)a  n  =  0,1,2,..,     , 

For  n  =  N,  simultaneous  periodic  solutions  will  occur  but  the  leading 

coefficients  of  the  trigonometric  series  will  vanish  up  to  and  including 

A^„  and  B^.,.  The  solutions  obtained  are 
^N        2N 

00 

(11)  u  ■  >  Ap  cos  2nx 

N+1  '^^ 

00 

(12)  u  =  r"  A_  sin  2nx  . 

For  the  other  alternative,  that  is 

(13)  c  «  -2na         n  -  0,1,2,... 

simultaneous  periodic  solutions  will  also  occur.  For  n=N,  only  the  leading 


Where  the  gamma-function  has  simple  poles.  See  [U6],  p.  236, 
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coefficients  up  to  A_,,  and  B„.,  will  differ.  The  solutions  obtained 

2N  ^N 


are 


N  00 

(lU)         u  -  T~  A     COS  2nx  +  T~     A     cos  2nx 


and 
(15) 


00 


u  "  r~  B     sin  2nx  +  P'     A 

r  "  fei  '^ 


sin  2nx 


The  values  of     b     for  the  solutions  with  period     n     are  obtained 
by  equating  the  complete     (N+l)-th     denominator  of    (7.9)  to  zero,   hence 

(le)         b  -  (2N+2)^  +   ...     . 


Similar  results  are  obtained  for  the  solutions  Co  and  So. 
Rewriting  the  recurrence  relationships  , 

(17)  [(2n*3)a-c]A2„^3/A2^^^  -[^^n-Da+cjA^^^^/A^^^^  -  2  [b  -  (2n.l)2j 

(18)  [(2n+3)a-c]B2„^3/B2„^l  "[(^n-l  )a+c]B2^.^/B2„^l  =  2  [b  -  (2n*l)2j 

As  before,    equate  the  left-hand  sides,   multiply  by     Ap     ,   ^2n+l* 
divide  by     2a     and  rearrange  terms   to  get 

■'2n+l  ■  ^2n+3^2n+lJ 


(19)  [(n-3/2).(c/2a)j[A2^^3B2„^^-B2„^3A2„,J      - 

-   C(n-l/2)-(c/2a)][A2^^^B2^_^  -  B^^^.A^^.,] 


or 


(20) 


*2n+3     ^2n+l 


^2n+3     ®2n+l 


(n-l/2)+(c/2a) 
(n+5>'2)-(c/2a) 


*2n+l     *2n-l 


^2n=«l     ^2n-l 


Continue  the  recursion  process  and  let 
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(20)     J^^^  -    (n-l/2)+(c/2a)     ;        I 


n+1 


(n*3/2)-(c/2a) 


then 


(21) 


*2n+3     ^2n+l 


^2n+3    ^2n+l 


2(.l)"  TT  CV^nVl 


and  as     n  ->  oo 


(22)  [P  (3/2  -  c/2a)/  p  (-1/2  +  c/2a)]A^B^  -  0  . 

Using  the  same  arguments  as  before,  either 

(23)  A.^B^  ■  0  or  (c/2a)-l/2  is  zero  or  a  negative  integer. 

If  A.^B-  -  0,  either  A.  «  0  and  all  the  A's  vanish  or  B 
and  all  the  B's  vanish.  Hence  simultaneous  periodic  solutions  do 
net  exist;  or 


(2h)  c  =  (2n+3)a, 


n  «  0,1,2,... 


and  simultaneous  periodic  solutions  do  occur.  The  leading  coefficients 
vanish  as  before,  and  for  n  ■  N-1  the  solutions  are 


00 

(2?)  u  -  r 

and 

00 


A2j^^j^  cos(2n+l)x 


A2j^^-j^  sin(2n+l)x     • 


(26)     u  -  p 

When  the  alternative  holds,  that  is 


(27)  c  -  -(2n-l)a  , 


n  "  0,1,2,.,, 


simultaneous  periodic  solutions  occur  with  their  leading  terms  differing. 

These  solutions  are  for  n  ■  N 

n-1  00 


(28)     u  »  ^  A^     .   cos(2n+l)x  +  T~  ^Pn+l  ^°s(2n+l)x 
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N-1  00 

(29)  u  -  5^  B-^^T   sin(2n+l)x  +  J^.  A2n+1  sin(2n+l)x  . 

o  N 

The  values  of  b  for  solutions  of  period  2n  are  obtained  by 
equating  the  complete  (N+l)-th  denominator  of  the  continued  fraction 
(7»18)  to  zeroj  hence 

(30)  b  -  (2N+1)^  +  ...  . 

As  a  convenience  for  future  reference,  the  group  u-solution 
pairs  (11)-(12),  (Ih)-(15),  (25)-(26),  and  (28)-(29),  will  be  referred 
to  as  Ince-type  solutions. 

These  results  may  now  be  summarized  in  the  followingt 

Theorem  8>1  -  The  Coexistence  Theoran  for  the  Whittaker  Equation. 

Two  periodic  solutions  of  period  n  or  2n  will  exist 

simultaneously  if  and  only  if  m  is  a  positive  or  negative  integer  or 

zero,  and  h  assumes  one  of  the  values  h   which  are  roots  of  the 
'  n 

characteristic  equation  for  h. 

Corollary  8«1 

If  m  is  an  odd  integer,  the  solutions  have  period  n  ;   and  if  m 
is  even,  the  solutions  have  period  2n. 

Theorem  8.2 

If  m  »  -  r,  r  =»  0, l,2,...j  and  h  "  h  ,  then  the  resulting 
Whittaker  equation  will  have  at  most  either 

(i)  r/2  odd  intervals  of  instability  for  m  even,  or 

(ii)  (r+l)/2  even  intervals  of  instability  for  m  odd. 

When  m  ■  0,  the  resulting  Whittaker  equation 

(31)   y"  +  (h  +  2q  cos  hx)j 
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is  a  special  case  of  Mathieu's  equation.  Theorem  8.2  implies  that  for 
m  »  0  all  the  odd  intervals  of  instability  vanish.  That  this  is  not  a 
contradiction*  may  be  seen  by  observing  that  the  least  period  of  this 
equation  (32)  is  ii/2. 

Theorem  (8.2)  indicates  that  for  r=l,  all  even  intervals  of  insta- 
bility, except  the  zero-th, vanish.  This  corresponds  to  m»l  and 

c»0  {  and  also  m»  -1  and  c-  -8q  (or  c-  -2a),  From  the  recurrence 

2 
relationship  we  see  that  b*0  or  h"  -2q   is  the  single  characteristic 

value  not  contained  in  the  h  •  For  all  other  values  of  b,  two 

solutions  of  period  n  coexist.  The  corresponding  Whittaker  equations 

are 


r   2  2  •^ 

(33)  y"  +  L~2q  +  Uq  cos  2x  +  2q  cos  Uxjy  -  0 

r-   2  2      •"! 

(3U)  y"  +  ].-2q  -  Uq  cos  2x  +  2q  cos  Uxjy  -  0 


This  vanishing  of  all  but  the  zero-th  even  interval  is  not  a  unique 
occurrence.  It  happens  again,  for  example,  in  the  case  of  the 
Frequency-modulation  equation. 

9.  The  impossibility  of  simultaneous  polynomial  solutions 

The  theorems  of  Section  7  indicate  the  possibility  of  solutions  with 
period  n  or  2n  in  which  the  trigonometric  series  terminate,  that  is 
to  say,  trigonometric  poljmomial  solutions  of  equation  (7«3)  may  occur. 
These  poljmomial  solutions  will  occur  if  and  only  if  c»  -na,  where  n 
is  any  positive  integer  0,1,2, ...etc.  The  coexistence  theorem  of  Section 
8  appears  to  imply  the  non-existence  of  simultaneous  polynomial  solutions, 


See  Borg's  result,  theorem  (2.1). 

WW 

The  polynomial  solutions  are  considered  in  Ince  ^9\< 
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since  the  values  of  b  requii^d  are  not  obtained  by  equating  one  of  the 
denominators  of  the  continued  fraction  to  zero. 

Consider  a  few  particular  examples  for  c-na  where  equation  (7.3) 
has  finite  solutions. 

If  0=0,  then  either  b-0  or  b  takes  the  values  given  by 
equation  (7.11).  If  b'O,  then  the  equation  to  be  solved  is 

(1)    u"  +  (a  sin  2x)u'  "  0 

and  the  solutions  are 


f 


(2)  u  ■  A  and        u  «  K 

^         0 


exp(:^  cos  2x)dx  * 


J 


The  second  solution  is,  of  course,  not  periodic. 

If  equation  (1)  is  transformed  to  the  normal  form  of  the  Whittaker 

equation,  the  degenerate  appearance  disappears  and  it  becomes 

2  2 

(3)    y"  +  (-2q  +  Uq  cos  2x  +  2q  cos  Ux)y  • 

If  c=  -2a,  then 

2         2 
(Ua)   u  ■  A  +  A^  cos  2x    for  b  -  Ub  -  Ua  »  0 

and  also 

(Ub)   u  -  Bp  sin  2x        for  b  ■  U  • 

If  c  ■  -Ua,  then 

(5a)   u  -  A  +-Ap  cos  2x  +  A.  cos  Ux 

for  b-  a/2)(2a-c)(c)/2^ 

1  -  b/2^  +  (t^a-c)(2a-»-c)/l6^ 


1  -bA^ 


and  also 
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(5b)    u  ■  B-  sin  2x  +  B,  sin  hx. 

1  -bA^ 


If  c  =  -a,  then 


(6a)    u  »  A^  cos  X     for  b  =  1  +  a 


and  also 
(6b)    u  -  B^  sin  x 

If  c  =  -  3a,  then 


for  b  -  1  -  a  , 


(7a)    u  =  A^  cos  X  +  A-  cos  3x 

for  b^  -  (10+2a)b  -  (3a^-l8a-9)  -  0 
and  also 

(7b)    u  ■  B^  sin  x  ♦  B^  sin  3x 

for  b^  -  (10-2a)b  -  (3a^+l8a-9)  »  0  . 

It  is  seen  that  for  a  given  value  of  c,  (a/O),  the  value  of  b 
is  different  for  the  even  and  odd  poljmomial  solutions.  It  is  clear, 
by  examining  the  continued  fractions  for  the  even  and  odd  solutions, 
that  the  same  pattern  will  continue  for  all  n. 

However,  the  impossibility  of  simultaneous  polynomial  solutions  can 
be  demonstrated  in  a  simple  and  direct  manner.  Assume  v  and  w  are 
two  solutions  of  equation  (7.3)  with  v(0)»l,  v«(0)»0,  w(0)=O,  and 
w«(0)-l,  then 

(8)  V*  +  (a  sin  x)v«  +  (b  +  c  cos  2x)v  -  0        and 

(9)  w"  +  (a  sin  x)w»  +  (b  +  c  cos  2x)w  ■  0, 


-  28  - 


Multiply  (8)  by  w  and  (9)  by  v,  and  subtract  (8)  from  (9),  then 

(10)  vw«  -  wv»  +  a  sin  2x  (vw»-wv»)  -  0  , 
Rewrite  as 

(11)  — i r—  »  -  a  sin  2x 

^  '       vw«  -  wv' 


and  integrate.  Hence 
log(vw»-wv')  -  ^  COS  2x  -  *• 


or 

V     ▼» 

(12) 

w    w' 

2 
■  exp(-  a  sin  x)  • 


If  V  and  w  are  both  finite  trigonometiric  polynomials,  then 
the  left-hand  side  of  equation  (12)  will  also  be  finite.  The  right-hand 
side,  for  a/O,  cannot  terminate,  and,  therefore,  v  and  w  cannot 
both  be  polynomial  solutions.  Therefore,  the  transformed  Whittaker 
equation  cannot  have  simultaneous  polynomial  solutions. 

The  question  of  whether  the  Whittaker  equation  can  have  coexisting 
polynomial  solutions  will  be  discussed  in  Chapter  V. 

10,  Imaginary  coefficients  and  coexistence 

If  the  coefficient  of  the  cos  Ux  term  in  the  Whittaker  equation 
is  taken  to  be  negative,  the  equation  is 

(1)  y"  +  (h  ♦  p  cos  2x  -  2q  cos  Ux)  . 
Assuming  solutions  of  the  form 

(2)  y  ■  u  exp(i  q  cos  2x)  , 

a  differential  equation  analogous  to  (7.3)  is  obtained 
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(3)    u*  +  (a  sin  2x)u»  +  (b  +  c  cos  2x)u  ■  0 

with 
(U)     a  -  -  Uiq,-    b  •  h-2q^;     c  «  p+Uiq  . 

From  the  results  obtained  in  Section  8,  it  can  be  concluded  that 
simultaneous  periodic  solutions  of  period  n  or  2n  can  occur  only  if 

(5)     p  "  hmiq 

where  in  is  any  integer  positive  or  negative. 

If  q  is  imaginary,  then  equation  (1)  returns  to  the  previous 
case.  However,  if  q  is  real,  p  must  be  imaginary  for  coexistence  of 
fundamental  periodic  solutions. 

For  Hill's  equation 

y«  +  [h  +  (2f(x)]y  -  0  , 

where  (^(x)  is  a  real  periodic  function  of  x  with  period  n,  it  is 
known  that  periodic  solutions  can  never  occur  if  h  is  complex  or 
imaginary  (see  ^32],  pp.  9-11).  This  result,  however,  does  not  exclude 
the  possibility  of  periodic  solutions  if  (2f(x)  is  complex. 

The  existence  of  periodic  solutions  of  (1)  can  be  inferred  from  (2) 
and  (3),  Equation  (2)  shows  that  the  solutions  of  (1)  are  periodic 
whenever  u  is  periodic »  From  Section  7  it  can  be  seen  that  the 
periodicity  of  the  solutions  of  equation  (3)  is  undisturbed  by  making  a 
and  c  imaginary,  although  the  coefficients  of  the  infinite  series, 
the  A's  and  B»3  may,  in  general,  be  expected  to  be  complex. 

From  the  tabulation  given  by  Klotter  and  Kotowski*  values  of  h 


See  [28],  Table  3,  p.l5U.  Note  the  factor  of  h  required  by  difference 
between  equation  (1)  and  their  equation.  The  corresponding  tabulation 
for  the  Whittaker  equation  (7.1)  is  given  in  their  Table  2, also  on  p.  15U. 


-  30  - 


corresponding  to  the  first  five  periodic  solutions  of  equation  (1)  can 
be  obtained  for  a  limited  range  of  2q  (0  to  2)  and  p  (0  to  6), 
The  values  of  h  are  real,  and  range  between  -2.8  and  6.3   The 
tables  give  only  positive  values  for  p.  However,  the  stability  chart 
is  symmetrical  around  p»0.  This  fact  can  be  seen  at  once  by  making 
the  substitution  x  «  -  Xo 

11,  An  extension  of  the  Whittaker  equation 

The  Whittaker  equation  has  three  arbitrary  parameters.  As  noted 
in  Section  5  equations  with  more  parameters  could  be  treated  if  they 
could  be  transformed  into  a  Whittaker  equation.  The  extra  parameters 
are  absorbed  in  the  transformation.  As  an  example  of  this,  consider 
the  extension  of  the  Whittaker  equation 

(1)  y"  ■•■  (hr   sin  2x)y'  +  (k  +  p  cos  2x  +  s  cos  Ux)y  =  0  . 

This  is  readily  transformed  into  the  standard  form  of  the  Whittaker 
equation  by 

(2)  y  =  w  exp(r  cos  2x) 
and  the  equation  in  w  is 

(3)  w«  +  [(k-2r  )  +  (p-Ur)  cos  2x  +(s+2r^)  cos  Ux]w  -  0  . 

Equation  (3)  is  transformed  into  equation  (7«3)  by  the  transforma- 
tion  (7.2)  with 

(U)  q-i   [(s  *  2r2)/2j2 

while  equation   (l)  is  transformed  into   (7.3)  by  (7.2)  with 
(5)  q  -  r  :^   f(s  +2r2)/2l2   « 
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An  equation  similar  to  (1)  which  also  reduces  to  the  Whittaker 
equation  is 

(6)  y"  +  (hr  cos  2x)y'  +  (h  +  p  cos  2x  -  Ur  sin  2x  +  q  cos  Ux)y  ■  0  . 

Equations  (l)  and  (6)  have  four  parameters,  and  combining  them 
leads  to  an  equation  vdth  five  parameters 

(7)  y"  +  (Ur  cos  2x  +  Us  sin  2x)y' 

+  (h  +  p  cos  2x  -  Ur  sin  2x  +  q  cos  Ux  +  Usr  sin  Ux)y  ■  0 
and  the  transformation  to  bring  equation  (7)  into  a  Whittaker  equation  is 

(8)  y  •  u  exp(s  cos  2x  -  r  sin  2x)  . 

Equation  (8),  with  8=0,  serves  for  equation  (6), 

2  2 
If  q»2(r  -s  )  in  equation  (7),  then  it  transforms  into  the 

2  2 
Mathieu  equation.  If,  in  addition  to  q"2(r  -s  ),  p"Us,  then  (7) 

transforms  into  the  equation  of  the  simple  oscillator 

(9)  y»  +  ky  -  0  . 

In  general,  if  P(x)  is  a  periodic  function,  it  is  always  possible 
to  generate  a  new  equation  by  setting 

(10)  y  =  u  P(x) 

in  the  Whittaker  equation.  Obviously  equations  so  generated  transform 
into  the  Whittaker  equation. 
The  equation 

(11)  y"  +  (h  +  p  sin  2x  +  q  cos  2x)y  -  0 

is  transfoimed  into  the  Whittaker  equation  by  a  simple  change  of 
independent  variable 
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(12)  X  "  t  -  n/U  . 
Similarly  the  equation 

(13)  y"  ♦  (h  +  p  cos  2x  +  p  sin  2x  +  q  sin  Ux)y  -  0 
is  transformed  into  the  Whittaker  equation  by 

(lU)    X  -  t  -  n/8  . 

Although  a  rnanber  of  simple  transformations  have  been  indicated 
here,  it  is  not  always  possible  to  find  a  transformation.  For  example, 
the  equation*" 

(15)    y"  +  (h  +  p  cos  2x  +  q  sin  Ux)y  -  0 

cannot  be  transformed  to  a  Whittaker  equation. 


See  [28],  p.  150.  The  least  period  has  been  changed  from  2n  to  n. 
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III,  The  Ince  Equation 

12,  Introductory  remarks 

An  equation  which  has  four  arbitrary  parameters  and  is  thus  more 
general  than  the  transformed  Whittaker  equation  (7.3),  but  which  also 
leads  directly  without  transformation  to  three-terra  recurrence 
relationships,  is 

(1)   (1  +  a  cos  2x)y»'  +  (b  sin  2x)y'  +  (c  +  d  cos  2x)y  -  0  , 

It  will  be  assumed  that   |a|  <  1,  setting  aside  the  question  of  what 
happens  when  |a|  >  1,  and  the  equation  has  singularities  on  the 
finite  portion  of  the  real  axis.  With  this  restriction,  (1)  is  closely 
connected  to  some  known  equations  with  elliptic  functions  as  their 
periodic  coefficients.  This  connection  will  be  explored  in  detail  in 
Chapter  IV,  As  an  example,  the  important  Lame  equation  can  be  trans- 
formed  into  a  special  case  of  equation  (1), 

The  Ince  equation  is  interesting  because  it  permits  a  discussion 
of  the  coexistence  problem  in  a  case  with  four  arbitrary  parameters, 
which  is  an  extension  of  previously  treated  cases.  It  is  not  difficult 
to  see  by  inspection  that  (1)  is  the  most  general  non-singular 
trigonometric  equation  which  yields  three-term  recurrence  formulas. 

With  certain  particular  values  for  the  parameters,  equation  (l) 
has  application  to  electric  circuit  theory.  For  example,  with  b"d»0, 
the  Ince  equation  becomes  the  well-known  Frequency-modulation  equation, 

13,  Conditions  for  coexistence  of  periodic  solutions 

Again  we  examine  the  question  of  simultaneous  existence  of 


With  b-  -a.  See  [2UJ,  p.  8U. 
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periodic  solutions  by  assioming  the  following  four  solutions 

00 


(Ce) 

y 

-  y~  A„     cos  2nx 

(Se) 

y 

'  f  ^2n  ^^"  2nx 

(Co) 

y 

00 

"EAjn+l  cos(2n^l)x 
o 

(So) 

y 

00 

■  JZ  ^2n+l  sin(2n+l)x  • 

0 

Consider  the  solution  Ce.  By  direct  substitution  in  equation 
(1),  the  recurrence  relationships  are: 

(|  -  b  -  2a)A2  +  cA^  -  0 

(1)  (|  -  2b  -  8a)Aj^  +  (c  -  U)A2  +  dA^  -  0 

(|  -  3b  -  l8a)A^  +  (c  -  16)A^  +   (|  +  b  -  2a)A2  -  0 
and  the  general  tem,  for    n  >  2,  is 

(2)  (d  -(2n+2)b  -  (2n+2)^a2A2j^  +2[c  -   (2n)^A2^ 

+  [d  +  (2n-2)b  -  (2n-2  )^anA2^_2  "  °  • 

Now  assume  a  solution  Se  also  exists.  Then  the  corresponding 
recurrence  relationships  are: 

(|  -  2b  -  8a)Bj^  +  (c  -  U)B2  -  0 

(3)  (|  -  3b  -  18a )B^  +  (c  -  16)B^  ♦  (|  +  b  -  2a)B2  -  0 
(|  -  Ub  -  32a)Bg  +  (c  -  36)B^  ♦  (|  +  2b  -  8a)Bj^  -  0 
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and  the  general  term,  for  n  >  2,  is 


(U)        [d  -   (2n+2)b  -   (2n+2)2aj\         +2[c  -   (2n)^jB 


2n 


+  [d  +   (2n-2)b  -  (2n-2)^a|B2^_2  "  ° 


Rewrite  equations    (2)  and   (U)  as 

A  A 

[(2n+2)^a  +  (2n+2)b  -  d]  4^    +  [(2n-2)^a  -  (2n-2)b-d]  -l^ 

*2n  ^2n 


2[c  -  (2nrj 


[(2n+2)^a  +  (2n+2)b  -  dj  4^  *  [(2n-2)^a  -  (2n-2)b  -  dj     ^""^ 


2n 


"2n 


-  2[c  -  (2n)2]    . 

Equating  the  left-hand  sides  of  these  equations,   clearing  fractions 
and  rearranging  terms,   we  get 


(5) 


A  A 

*2n+2     '^2n 


^2n+2     ^2n 


f2n-2)  a  -  (2n-2)b  -  d 
(2n+2)  a  +  (2n+2)b  -d 


^2n     ^2n-2 


^2n    ^2n-2 


Adopt  the  notation 


(6)         K„  -  (n-X)2  -  (n-l)  |j  -  gj 


and 


>>n.l  ■  (n*l)'*(n*l)  h  "  ^ 


whence 


(7) 


B 


n+2     ^2n 


2n+2       2a 


n 


2  TT     (K^^r^A^Bg     • 
r-1 
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Since  unlimited  differentiability  of  the  trigonometric  series  requires  that, 
for  any  K  >  0, 


(3)      Lim 

n->oo 


^n+2  *2n 


^2n+2  ^2n 


n  ■  0 


equation  (8)  implies  that 

(9)  I  (1  ^  p  ^  q)  I  (1  ^  p  -  q)  A  B  -  0 
I  (-  P  ♦  q)   I  (  -  P  -  q)   °  ^ 

where  p  -  bAa,  Q  -  (b^  ♦  Uad)-^/^Aa  . 

The  coexistence  question  is  settled  for  the  solutions  of  period  n    by  an 

examination  of  equation  (9)»  Briefly  stated,  unless  the  argument  of  one  of  the 

two  gamma-functions  in  the  denominator  of  (9)  is  zero  or  a  negative  integer,  two 

solutions  of  period  n     cannot  exist  simultaneously  unless  A  Bp  ■  0.  In  this  case, 

the  recurrence  relations  (2)  and  (U)  show,  that 

(2n+2)^a  +  (2n+2)b  -  d  »  0 

for  at  least  one  integer  n  »  1,2,3,  ...  t   unless  either  all  A^  ■  0  or  all  B-  "0. 

This  is  equivalent  to  stating  that  at  least  one  of  the  Gamma  factors  in  the 

numerator  of  (9)  must  be  infinite.  Altogether,  we  have  as  a  necessary  condition 

for  coexistence  of  two  linearly  independent  solutions  of  period  n     that 

(10)  p  +  q-in,     n"  0,1,2,... 

The  characteristic  values  for  c  sure  taken  from  the  continued  fraction 
obtained  in  the  usual  way.  In  this  case,  however,  there  are  three  possibilities. 
The  values  of  d  may  be  taken  from  a  terminated  continued  fraction  or  from 
equating  the  appropriate  infinite  complete  denominator  to  zero,  as  in  the  Whittaker 
case.  In  addition,  the  values  of  d  may  be  taken  from  the  terminated  continued 
fraction  which  begins  with  some  nth  denominator. 

Now  consider  the  solution  C  .  Again,  by  direct  substitution,  the  recurrence 
telationshipe  are 
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(d-3b-9a)A-   ♦   [(2c-2)+d-b-a]A^  -  0 

(11)  (d-?b-25a)Aw  +  2(c-9)A^  +   (d*b-a)A^  -  0 
(d-7b-U9a)A^  ♦  2(c-25)A^  ♦  (d+3b-9a)A^  -  0 

and  the  general  term  for    n  >  1,  is 

(12)  Id  -(2n+3)b  -(2n+3)^a]A2^^^  ♦  2[c  -(2n+l)^]A2j^^^ 

♦  [d  +(2n-l)b  -(2n-l)^a]A2j^_^  -  0  . 
Similarly,  for  the  solution  So,  the  recurrence  relationships  are 
(d-3b-9a)B^  ♦  [2(c-l)-d+b+a"jB.  -  0 

(13)  (d-5b-25a)B^  +  2(0-9)33  ♦  (d+b-a)B^  -  0 
(d-7b-U9a)B  ♦  2(c-25)B^  +  (d+3b-9a)B3  -  0 

and  the  general  term  for  n  >  1,  is 
(lU)   [d  -(2n+3)b  -(2n+3)^a]B2^^3  +  2[c  -(2n+l)^]B2^^^ 

♦  [d  ♦(2n-l)b  -(2n-l)^ajB2^^^  -  0  . 
Rewriting  equations  (12)  and  (ik)   as 

I(2n+3)^a  +(2n+3)b-d]  J^  ♦  [(2n-l)^a  -(2n-l)b-dJ  J^  -  2[c-(2n+l)^]  , 

*2n*l  *2n+l 

r(2n>3)2a  *(2n*3)b-d]  l^  *   r(2n-l)2a  -(2n.l)b-d]  ^  -  2Ec-(2n*l)2j  . 
•-  «2n+l  2a+l 
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Equating  the  left-hand  sides  of  these  equations,  clearing  fractions  and 
rearranging  terms,  we  get 


(15) 


^n+3  ^2n+3 


*2n+l  ^2n+l 


(2n-l)_a.(2n-l)b-d 
(2n+3)^a+(2n+3)b-d 


*2n*l  ^2n+l 


2n-l   2n-l 


Adopt  the  notation 


(16)        J^^^  -  (2n-l)^-(2n-l)  |  -  |      and 


n+1 


(2n+3)^+(2n+3)  |-| 


whence 


(17) 


A  B 

2n+3      2n+3 


*2n+l    ®2n*l 


n+1 

"  2ir  (v^n>vi  • 


Again,   v:e  'mow  that   f;>r  h.t/  k  >  0 


(18) 


Lim 
n->oo 


*2n+3     ®2n+3 


*2n+l    ^2n+l 


^  ^  n 
n     •  0 


and,   therefore, 

(19)       n3/2  -^  P  ^  q)  r  (3/2  ^  p  ■  q^        A,B,  -  0 

n-  1/2  -  p  ♦  q)   r  (-  1/2  .  p  -  q)       ^  ^ 

2  X/2 

where  again    p  ■  bA^j     q  °  (b     +  had)  '     Aa     • 

The  coexistence  question  is  settled  for  the  solutions  of  period  2n 

by  an  examination  of  equation  (19).  Again,  unless  the  argument  of  one  of  the 

four  gamma-functions  in  (19)  is  zero  or  a  negative  integer,  two  solutions  of 

period  2n  cannot  exist  simultaneously;  or  what  is  equivalent 


(20)   ptq"in-l/2;     n-  0,1,2,...  • 
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Here  again  there  are  three  vrays  in  which  the  values  of  c  may  be 
obtained  from  the  continued  fraction,  as  stated  in  the  remark  following 
equation  (10),  It  may  be  noted  that  for  the  values  of  c  obtained  from 
the  terminated  continued  fraction  beginning  with  the  initial  tenns, 
coexistence  does  not  occur,  but  can  occur  in  the  other  two  instanceso 

Theorem  13 ♦  3--The  Coexistence  Theorem  for  the  Ince  Equation 

Simultaneous  periodic  solutions  will  exist  for  the  Ince  equation  if 

and  only  if*: 
(i)   For  solutions  with  period  n, 

P  *  ^  =  i  \  f  "i*'^  "  0>1>2,  ... 

or 


p  -  q  -  +  n^ 
(ii)   For  solutions  with  period  2n 

p  +  q  -  ±  n^  -  1/2  ;    n^jn^  -  0,1,2,  ... 
or 

P  -  q  -  ±  n2  -  1/2  . 

Corrolary  13.1 

There  are,  in  general,  three  types  of  simultaneous  periodic  solutions, 
(The  m  used  below  stands  for  either  (2n)  or  (2n+l)  accordingly  as 
the  solutions  are  of  period  n  or  2n)» 

(i)    If  in  theorem  (13.1)  n  <  0,  then  the  solutions  have  the  form 
(tiruncated  solutions) 

00 

y  ■  y*"  A  COS  mx 


■I^*m 


00 

y  "  r~  A  sin  mx 


It  is  assumed  that  c  has  the  required  characteristic  value. 
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(ii)   and  if  n  ^0»  then  the  solutions  have  the  form  (split  solutions) 
N-1  oo 


y  -  ^  A^  cos  mx  ♦  ^  A^ 

N-1  00 

y  ■  y  B  sin  mx  ♦  T"  A 


(iii)     and  finally  the  solutions  which  result  irtien  more  than  one  of  the 
gajnma-functions  in  equations   (9)  and  (19)  has  a  simple  pole.     If,  for  the 
same  values  of  the  parametei^  of  the  equation, 

p  -  q  ■  -r-1    and  also     p  -  q  ■  s       or 

p  -  q  »  -r-1/2   and  also     p  _  q  »  s  -  1/2 

where  r  and  s  are  positive  integers,  then  simultaneous  polyncanial 
solutions  will  occvir  ^ich  may  be 

(a)  truncated  for  r  <  s 

M 

y  ■  5~  A  cos  mx      0  <  N,M  <  oo 

r  ° 

M 

y  ■  r~  A  sin  mx 

r  " 

(b)  split  for  r  >  s 

n-1  M 

y  ■  T~  A  cos  mx  +  7""  A  cos  mx     0  <  N,K,  <  oo 

n-1  M 


sin  mx  • 


Corollary  13*2 

The  Ince  equation  will  have  at  most 
(i)  (n+1)  even  intervals  of  instability  if 

piq-*nj     n-  0,1,2,  ... 


-  Ul  - 


(ii)       I  p.  I     odd  intervals  of  instability  if 

p+q-in-l/2  n+  0,1,2,    ... 

(iii)  (2ri+l)  even  and  odd  intervals  of  instability  (i.e.,  only  a  finite 
number  of  intervals  of  instability  remain)  if  simultaneously 
P±q-±n,   p*q-±n-  1/2,   n  -  0,1,2,  ...  . 

Corollaiy  13.3 

(i)  For  d  "■  0,  all  even  intervals  of  instability,  except  the  zero-th 
interval,  vanish. 

(ii)   For  d  "  a  +  b,  all  odd  intervals  of  instability  vanish. 

(iii)  For  d  =  0,  and  a  ■  -b,  all  intervals  of  instability  except  the 
zero-th,  vanish. 

The  discussion  of  (ii)  and  (iii)  as  related  to  Borg's  theorems  (2.1) 
and  (2.2)  is  deferred  until  Chapter  IV. 

Ih.   Special  cases.  The  Frequency-modulation  equation 

An  interesting  special  case  of  the  Ince  equation  is  the  Frequency- 
TTodulation  equation  which  has  important  applications  in  electronic 
communications.*  It  is  the  equation  for  a  sertes  inductance-capacitance 
electric  circuit  with  no  resistance  in  which  the  capacitance  varies 
sinusoidally,  and  has  been  discussed  by  McLachlan^r  J  and  treated  extensively 
by  Cambi  ^ -"  *  L^ '  .  The  equation  is  written 

(1)   (1  ♦  a  cos  2x)y"  ♦  cy  -  0 

*nd  again  the  restriction  that  |a|  <  1  is  imposed  to  avoid  singularities. 
Cambi  considers  the  cases  where  |a|  >  1,  but  a  discussion  of  this  will  be 
deferred. 


The  literature  is  voluminous.  See,  for  example,  [^  and  [Ut] . 
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Since  equation  (1)  is  the  Ince  equation  vjith  b  ■  0,  and  d  •  0, 
the  results  for  the  Ince  equation  are  directly  applicable.  Using  these 
values  for  b  and  d  we  obtain 

Theorem  m«3 

In  the  Frequency-modulation  equation,   all  the  even  intervals  of 
instability,   except  the  zero-th,   always  vanish,   and  none  of  the  odd 
intervals  of  instability  ever  vanish* 

The  theorem  follows  at  once  from   (13.7), 

'^2n+2  ^2n  '  ^2n  ^2n+2  *  ° 
and  (13.3  7)       A^B^  -  0     • 

Two  additional  interesting  equations  can  be  evolved  from  equation  (1)« 
They  are 

(2)  (1  ♦  a  cos  2x)v"  -  (2a  sin  2x)v'  +  cv  -  0 
ana 

(3)  (1  +  a  cos  2x)w«  -  (ha  sin  2x)w»  ♦  (c  -  Ua  cos  2x)w  ■»  0 
and  are  characterized  by  the  fact  that  if 

y  «  F(x,r),  then  v  ■  F'(x,a),  and  w  -  F"  (x,a)  , 

From   the  recurrence  relationships  for  the  Ince  equation,  it  can  be  seen 
that  theorem  (itt.l)  also  holds  for  equations  (2)  and  (3).  The  first 
recurrence  relation  for  the  Ce  solution  shows  that  the  constant  term  A 
which  is  present  for  (1)  is  zero  for  the  derived  equations  as  expected. 

Equation  (2)  may  be  considered  to  be  the  equation  for  a 
resistanceless  series  inductance-capacitance  circuit  in  which  the  inductance 
is  varying  periodically. 
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A  more  general  form,  of  which  equation  (2)  is  a  particular  case, is 

(U)  (1  ♦  a  cos  2x)y"  ♦  (b  sin  2x)y'  ♦  cy  «  0  . 

Here  again,   all  of  the  even  intervals  of  instability,   except  the  zero"th 
interval,   always  vanish.     But  the  odd  intervals  can  also  vanish,    since  for 

(5a)     b  «  (2n-l)a  or 

(^)     b  «  -  (2n+3)a 

simultaneous  solutions  with  period  2n  can  occur.  Equation  (U)  is 
perhaps  the  simplest  in  which  it  can  happen  that  only  a  finite  number  of 
intervals  of  instability  remain. 

The  periodic  solutions  of  equation  (I4),  aside  frcm  the  fact  that 
those  with  period  n  occur  almost  always  in  even  and  odd  pairs,  behave  in 
a  way  similar  to  the  solutions  of  the  VjMttaker  equation.  The  solutions  with 
period  n  can  terminate  if 

(6a)     b  -  2(n  -  l)a 

and  finite  polyncanial  solutions  can  occur.  If 

(6b)     b  »  -  2(n+l)a 

the  coefficients  of  the  trigonometric  series  up  to  and  including  A_   and 
Bp  can  vanish  identically. 

In  the  same  manner  that  equation  (3)  was  obtained  from  equation  (2), 
equation  (U)  leads  to 

(7)   (1  +  a  cos  2x)u»  +  ([(b-2a)  sin  2x]u»  ♦  (c  +  2b  cos  2x)u  ■  0 

which  is  a  particular  case  of  the  more  general  Ince  equation.  Again,  if  in 
(U) 

y  -  G(x),    then   u  -  0'(x) 


'hh' 


and  the  properties  of  its  solutions  can  be  expected  to  be  similar  to  those  of 
equation  (U)» 

If  we  set  b«2a,  then  (7)  reduces  to 

(8)  (1  ♦  a  cos  2x)u"  +  (c  +  Ua  cos  2x)u  »  0 
which  is  a  particular  case  of 

(9)  (1  ♦  a  cos  2x)y«»  +  (c  +  d  cos  2x)y  «  0  , 

This  last  equation  is  interesting  because  the  criterion  for  coexistence 
is  simply 

(10)  d  /-  n^a 

and  also  because  in  this  case  the  coefficients  of  A   and  B   vanish  in 

n       n 

two  of  the  three  recurrence  relations  in  which  they  appear.  The  usual 
alternative  conditions  (e.g.,  b  ■  2(n-l)a  and  b  +  -2(n+l)a,  etc.) 
occur  ,  so  to  speak^ simultaneously. 

Even  and  odd  intervals  of  instability  do  not  vanish  at  the  same  time, 
since  for  n  even,  even  intervals  of  instability  vanish;  and  for  n  odd, 
odd  intervals  of  instability  vanish. 

The  Frequency-modulation  equation  is  the  particular  case  of  equation 

2  2 

(9)  for  n  "  0.  If  we  set  c  »  n   whenever  d  «  n  a  ,  then  the  equation 

(9)  reduces  to  the  equation  for  the  simple  oscillator,  and  the  solutions 

are  always 

A  cos  nx    and   B  sin  nx  . 
n  n 

An  indication  of  a  physical  application  for  equation  (9)  is  given 

in  the  next  paragraph  where  its  connection  with  equation  (15»7)  is 

described. 
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3  5^.   Normal  forms 


Rewrite  the  Ince  equation  as 


(1) 


b  sin  2x 


y«     ♦    ■"  —  —     .       yt  4   f^/a  ♦  c  -   (d/a)        1 
■^  1  +  a  cos  2x    ^         [_'         r+acosTxJ-' 


0  . 


Make  the  trcmsformation 

sin  2x 


(2)   y  -  u  exp(-  1/2 
and  obtain 


f  b  si 

J  ^ 


a  cos  2x 


dx) 


(3)   u«  ♦ 


A/.  ^  ^  '   d/a  -  b    (b^+Uab)sin^  2x 
^/^  *   1  ^  a  cos  Sx  "  TTT" TT? 


U(l  ♦  a  cos  2x)' 


u  •  0 


which,  with  a  little  manipulation,  becwnes 

ih)       (1  •»■  a  cos  2x)  u"  +  (h  +  p  cos  2x  +  q  cos  Ux)u  "  0 
where  in  terms  of  the  parameters  of  the  Ince  equation 


h  -  (d/2)  -  (b'^/S)  -  (ab/2)  ♦  c  -  b 

p  «  d  ♦  (c  -  b)a 

q  -  (ad/2)  ♦  (b^/8)  +  (ab/2)  • 


Equation  (U)  is  the  so-called  normal  form  of  the  Ince  equation,  and  the 
results  for  the  Ince  equation  can  be  applied. 
Recalling  that 

equation  (2)  may  be  expressed  as 
(?)  u  -  yd  +  a  cos  2x)"  ^/^*  . 

Of  the  equations  in  Section  lU,  the  Frequency-modulation  is  in 
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normal  form,  and  the  normal  forms  for  the  others  appear  different  only  in 
the  absence  of  a  parameter.  For  equation  (lU.U),  where  d  »  0,  the  normal 
form  is 

(6)  (1  +  a  cos  2x)  u"  +  [(r-s)  ♦  ra  cos  2x  •♦•  s  cos  Uxjti  ■  0 

2 
with  r  "  c  -  b  and   s  »  b  /8  ♦  ab/2 

and  the  normal  form  for  equation  (lU.?),  where  b  ■  0,  is 

(1  +  a  cos  2x)  u"  +  [(d/2  +  c)  +  (d+ca)cos  2x  +  (ad/2)cos  Uxju  ■  0  • 
Closely  related  to  the  Frequency-modulation  equation  is  the  equation 

(7)  y"  ♦  2Ry«  +  [S/(1  +  a  cos  2x)]y  »  0 

which  describes  a  series  electric  circuit,  with  inductance,  capacitance, 
and  resistance,  in  which  the  capacitance  is  varying  periodically.  Applying 
the  transformation 

(8)  y  -  n  exp  [  (-R)dx  , 
equation  (?)  becomes 

(9)  (1  +  a  cos  2x)u»  +  [(S-R^)  -  aR^  cos  2xju  «  0 

which  is  equivalent  to  equation  (ltt«9)  with  a  modification  of  the  coefficients, 

16,   Oeneraliaations 

A  large  class  of  equations  with  periodic  coefficients,  which  can  be 
transforroed  into  the  Ince  equation,  and  whose  solutions,  consequently, 
behave  in  a  similar  manner,  may  be  generated  by  the  simple  expedient  of 
setting 
(1)   y  -  u  P(x) 

where  P(x)  is  a  periodic  function.  The  resulting  generalised  equation  is 
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(2)  (1  ♦  a  cos  2x)Pu"  +  [(b  sin  2x)P  +  2(1  +  a  cos  2x)P»]u' 

♦  [(c  +  d  cos  2x)P  ♦  (b  sin  2x)P«  +  (1  +  a  cos  2x)P«Ju  -  0  . 

It  is  relatively  easy  to  obtain  equations  of  the  form  in  (2)  from  the 
Ince  equation.  The  difficulty  lies  in  the  inverse  problem,  that  is, in 
recognizing  when  an  equation  is  in  the  form  of  (2),  and  can  therefore  be 
transformed  into  an  Ince  equation.  For  example,  if  in 

(3)  (1  +  a  cos  2x)u"  *  (s  sin  hx)u'  +  (h  ♦  p  cos  2x  ♦  q  cos  Ux  ♦  r  cos  6x)u  -  0 
the  parameters  have  the  values 

s  =»  -  ab 

h  •  c  -  b^/8  -  ab/2 

(U)         p  -  d  -  b  +  ab^A6 

q  -  b^/8  -  ab/2 

r  "  -  ab^A6  , 
then  the  transformation 

(5)  u  »  y  exp(b/2  cos  2x) 

reduces  (3)  to  an  Ince  equation. 

The  scheme  of  (1)  may  be  employed  to  introduce  additional  parameters. 
For  a  transformation  similar  to  ($)  of  the  form 

(6)  y  ■  u  exp(-  m  cos  2x) 

we  get  r-  .'■''  ■  ■'" 

(7)  (1  +  a  cos  2x)y"  ♦  [(b-2m)sin  2x  -  2am  sin  )ax]u» 

+  [h  +  p  cos  2x  ♦  q  cos  Ux  +  r  cos  6xju  "0 

with 
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2 
h  ■  c  -  inb/2  -  aa  +  a  /2 

2 
p»d-2in-t-am/U 

2 
q  ■  mb/2  -  ma  -  a  /2 

2 
r  »  am  A     • 

An  extension  of  the  Ince  equation  which  is  useful  as  well  as 
interesting  is 

(8)    (1  +  a  cos  2x)y"  +  (b  sin  2x)y'  +  (  c  +  d  cos  2x  +  ■-■ ^—^)j   -  o 

\  X  *  »  COS  cX  * 

Ccwiparing  this  with  equation  (l5,l),  it  is  seen  that  if  equation 

(8)  is  put  into  the  same  fortn  as  (l5»l),  it  will  differ  only  by  a  term 

2 
r/(l  +  a  cos  2x)  •  It  is  therefore  natural  to  consider  using  the 

transformation  (15.2).  The  resulting  equation  is 

2.,..x_.    2 


u"    ••■   jd/a  + 


(9)  u»    *  Id/a  *  r^/^-^       ■     ♦  Ur-(b  ^Uab)sin     2x 

^  ^  '     ^  ♦  a  cos  2x      TTrrm  ^rr? 


u 


U(l  ♦  a  cos  2x) _ 

which  is  alaost  identical  with  equation  (1^.3 ),  ''nd  can  easily  be  wrttten 
in  the  same  form  as  equation  (15. U),  the  normal  form  of  the  Ince  equation, 

(10)  (1  +  a  cos  2x)  u»  ♦  [h  +  p  cos  2x  +  q  cos  Ux]u  »  0  . 

This  same  transformation  (15«2)  will  reduce  equation  (8)  to  the 
Frequency-modulation  equation  if  the  parameters  in  (8)  are 

c  -  H  +  (b/2a)^ 

(11)  d  •  b  -  (b^/8a) 

r  -  (b/2a)^  +  (m^/8)  . 

Transformations  of  the  Ince  equation  emplojring  elU  ptic  functions 
are  also  interesting,  and  lead  to  known  equations.  Some  of  these 
transformations  vdll  be  discussed  in  Chapter  IV, 
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IV.  The  Herndte  Elliptic  Equation 
17 •  Intyodpctoyy  remnrks 

The  theory  of  differential  equations  with  doubly-peidodic  coefficients, 
developed  by  Hemite  and  Picard,  preceded  the  Floquet  theory  of  differential 
equations  with  singly-periodic  coefficients.  Perhaps  the  best  known 
equation  witlv  elliptic  functions  as  coefficients  is  the  Lame'  equation,  for 
which  there  is  an  extensive  literature.  Nonetheless,  the  coexistence 
question  for  the  Lame'  equation  was  not  discussed  until  19U0  when  Ince 
transformed  it  into  a  trigonometric  equation.  It  is  interesting  to  note 
that  Ince  attributes  the  idea  for  his  paper  to  a  question  in  a  letter  of 

ri2*] 

Herroite'-  -■ .  The  idea  is  essentially  to  consider  sin  n(ara  z)  instead  of 
sn  nz  as  the  true  analogue  of  sin  nx.  As  Ince  points  out,  the  coexistence 
question  could  not  be  treated  using  the  series  of  powers  of  elliptic 
functions. 

The  Hemite  elliptic  equation  was  introduced  in  November,  187 9^ in  a 

letter  to  E.  Heine,  and  a  slightly  simpler  equation  of  the  same  type  was 

flfl 
treated  earlier  by  Picard''  -' ,     The  normal  form  of  this  equation  was  called 

[20] 
the  Associated  Lame  equation  by  Ince  .  In  his  letter  to  Heine,  Hermite 

also  introduces  two  other  equations  with  elliptic  fimctions  as  coefficients, 

the  three  equations  differing  only  in  the  coefficient  of  the  first 

derivative  term.  All  three  of  these  equations  are  special  cases  of  an 

equation  introduced  bv  Darboux,  but  usually  referred  to  a^  Heun's  equation.  (see[7j 
and  [13],  PP»  161-196.) 

See  (211],  pp.  83-99.  Ince  refers  to  a  letter  of  Hermite  to  Dini,  Oeuvres, 

V  U,  p.  lUU. 
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Extensions  in  two  directions  are  indicated.  The  first  extends  the 
Hermite  elliptic  equation  in  a  manner  analogous  to  (16.8).  This 
equation  has  not  been  identified  with  a  previously  discussed  equation. 
The  second  extends  the  Lame'  equation  to  the  Lame  wave  equation  or 
Ellipsoidal  wave  equation,  and  to  the  Ultra-ellipsoidal  wave  equation. 
The  Ince  equation  remains  the  basis  for  these  discussions. 

The  point  of  departure  in  this  chapter  is  different  from  those  in 
the  preceding  ones.  In  the  case  of  the  Whittaker  equation  and  the 
extensions  of  the  Ince  equation,  the  basic  tool  was  the  transformation  of 
the  dependent  variable.  In  this  chapter  we  anploy  a  transformation  of 
independent  variable  to  change  from  elliptic  functions  to  trigonometric 
functions, 

18.   The  Hermite  elliptic  equation 

The  Hermite  elliptic  equation  which  contains  four  arbitrary  parameters 
is  written 

(1)  y.  *  iCi-ljk^s^j^cn^  y,  ^  0,.(„.r)(m.r*l)k2  sn^  xjy  -  0  . 

The  transformation  to  the  Ince  equation  is  accomplished  by  letting 

(2)  X  -  am(t,k) 

(3)  dt/dx  -  (1  -  k^  sin^  t)^^^ 

2 
where,  of  course,  k  <  1. 

Since  y»  -  yt«    and  y"  -  y(t«)^  +  yt« 

where  the  dots  indicate  differentiation  with  respect  to  t,  the  equation 

becojnes 

(U)  (1-k^  sin^  t)y  +  f(2r+l)k^  sin  t  cos  t]y  +  [h  -(m-r)(m+r+l)k^  sin^  t]y  -  0  . 


with 


(?) 
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Equation  (U)  is  the  Ince  equation 
(1  ♦  a  cos  2t)y  +  (b  sin  2y)y  +  (c  ♦  d  cos  2t)y  -  0 


a  -  k^/(2-k^) 

b  "  (2r+l)k^/(2-k^) 


c  =  [2h  -(ai-r)(m+r+l)k^]/(2-k^) 
d  -  (m-r)(m+r+l)k^/(2-k^) 

and  we  note  that  k  <  1  implies  (a(  <  1  which  was  the  restriction 
imposed  in  Section  12. 

With  these  values  for  the  parameters,  the  p  and  q  of  Section  13 
become 
(6)   p  =  (2r+l)A    and   q  -  (2m+l)/U 

Theorem  18 «] 

The  condition  for  coexistence  of  periodic  solutions  for  the  Hermite 
elliptic  equation  is  that  (r-nn)  or  (r-«)  be  an  integer.  If  (r+in) 
is  odd,  or  (r-n)  even,  the  solutions  will  have  period  2K,  while  if 
(r-Hn)  is  even,  or  (rnji)  odd,  the  solutions  will  hs-v^e  period  LK. 
The  Hermite  elliptic  equation  will  have  only  a  finite  number  of  intervals 
of  instability  if  (r-m)     and  (r^a)  are  simultaneously  odd  or  simultaneously 
even,  and  also  if  either  (r-*fli)  or  (r-m)  is  both  even  and  odd. 

If  m  and  r  are  integers,  then  (m+r)  and  (m-r)  are  odd  or 
even  together,  and  in  this  event  the  Hermite  elliptic  equation  will  have 
only  a  finite  nimiber  of  intervals  of  instability.  In  a  similar  way  the 
other  results  for  the  Ince  equation  may  be  translated  in  terms  of  r  and 
m. 


Corresponding  to  n  and  2n  in  the  trigonometric  form  of  the  equation. 
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The  explanation  of  Corallary  13.3  (lii)  is  now  easy.  The  criteria 
d  ■  0,  a  ■  -b  correspond  to  m  ■  r  and  r  «  -1,  and  therefore  (1) 
becomes  the  equation  for  the  simple  oscillator 

y  +  hy  -  0, 
thus  satisfying  the  conditions  of  Borg's  theorem  that  0(x)  must  be 
a  constant  for  all  intervals  of  instability  except  the  zero-th  to  vanish. 

If  r  «  m,  equation  (1)  becomes  the  Picard  elliptic  equation 

(7)  yn  ♦  '^^'^'^1,^;  ^  °^  ^  y.  *  hy  .  0 

which  is  transformed  by  (2)  into  the  special  case  of  the  Ince  equation 
given  in  equation  (lii.U)  with 

a  -  k^/(2-k^) 

(8)  b  -  (2r+l)k^/(2-k^) 
c  -  2h/(2-lc^)  . 

Therefore  p  -  q  =  0  and  consequently: 

Theorem  18.2 

For  the  Picard  elliptic  equation,  all  the  even  intervals  of  instability, 
except  the  zero-th  interval,  vanish  identically.     Odd  intervals  of 
instability  will  vanish  if,  and  only  if,     r    is  an  interger. 

The  Frequency-modulation  equation   (lii.l)  has  an  elliptic  fimction 
analog  in  the  special  case  of  the  Picard  elliptic  equation 

2 
/rt\  ••     .k    snxcnx_,   ^.      _- 

(9)  y"    ♦  ajf^j y   +  hy  -  0 


irtiere 
(10) 


a  -  b  -  k2/(2-k^) 
c  -  h/(2-k^)  • 
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Also 
(11)     r  -  0,   m  -  0 
vrtiLch  again  leads  to  the  results  of  theorem  Ih.l. 

The  normal  form  of  the  Hennite  elliptic  equation  is  obtained  by 
removing  the  term  in  y'  by  the  transformation 


(r*l)k  sn  X  en  X 


(12)   y.uexp[.J  — ^  ^ 
or    y  ■  u(dn  x)* 


dx 


] 


which  is  the  analog  of  equation  (15»2).  The  equation  obtained,  analogous 
to  (I5.h),  is 

(13)   y"  +  [h  ♦  p  sn  X  *   q(cn  x/dn  x)]y  »  0 

which  Ince  has  called  the  Associated  Lame  equation. 

An  extension  of  the  Hennite  elliptic  equation  to  an  equation  with 
five  parameters  is 


2 

/, ,  X    _  ^  2(r+l)k  sn  t  en  t   ,  ^ 
(lU)   y"  ♦  -i i-rrr 7'  + 


HiTT 


h  ♦  p  sin  X  +  -%■ 


dn  t 


y  ■  0 


which  transforms  directly  to  equation  (16,8). 

Equation  (13)  is  a  special  case  of  the  equation 


(15)   y" 


w  ^  /■  ^-1  ^l-2   2  ^  A   ^  B  dn^x  ^  C  k^  cn^:  .„ 
-h  ♦  m(m+l;k  sn  x  +  — 5-  +  5 —  ♦ 5 y 


T 


sn  x    en  X 


dn  X 


4 


which  DarboTUC  transformed  into  what  is  now  called  the  Heun  equation. 
An  interesting  equation  is 

(16)   y«  -  K  ^g^^^°"  ^  y  ♦  [a  ♦  (b  -  ak^)2n\   +  (c  -  bk^jsn^o  -  ck^  sn^tly  -  0 


which  transforms  into  the  Whittaker  equation. 
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(17)     y  +   (a  +  b  sin  t   +  c  sin  t)y  -  0, 

For     c  ■=  0,   equation   (ifc)  is  the  analog  of  the  Mathieu  equation, 

19,       The  Lsjne  equation 

The  Hemdte  elliptic  equation  reduces  to  the  Lame  equation  v^en 
r  •  -1 

(1)  y"     ♦  [h  -  m(m+l)k^  sn^xjy  -  0  , 

This  equation  contains  only  three  arbitrary  parameters.  The 
corresponding  trigonometric  equation  is 

(2)  (1  -  k^  sin\)y  -   (k^  sin  t  cos  t)y  +  |h  -  m(m+l)k^  sin^x|y  -  0 
which  is  the  Ince  equation  with 

a  -  -b  -  k^/(2-k^) 

(3)  c  -  i;2h  -  m(m*l)k^]/(2-k^) 
d  »  ra(m+l)k^/(2-k^)   o 

The  values  of     p     and     q     are 

(U)         p  -  -  lA  ,  q  -    (2m+l)A  • 

From     p+q  ■  n     or     p-q  -  n-]/2  ,   we  get     m  ■  2n  ,   and  from     p+q  ■  n-l/2 
and     p-q  ■  n  ,   we  get     ra  "  2n-l,     This  provides 

Theorem  19,1 

If  and  only  if  m  is  an  integer,  including  zero,  the  Lame  equation 
will  have  at  most  m+1  intervals  of  instability. 

Corollary  19,1 

If  and  only  if  m  is  an  interger,  then  solutions  of  the  Lame  equation 
with  period  2K  will  coexist,  and  those  wLth  period  UK  will  coexist 
for  the  same  value  of  m. 
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Since  the  Lame  equation  is  doubly-periodic,  it  has  fundamental 
imaginary  periods  in  addition  to  its  real  periods*  The  question  of 
coexistence  for  the  imaginary  periods  2iK«  and  UiK'  can  be  settled  by 
making  an  imaginary  transformation  of  the  equation  (seeQBjj and  [9], 
esp.  pp.  69-72).  Specifically, let 

(5)  X  -  i(x  -  K  -  iK«)  ,        H  -  m(m+l)  -  h 

and  the  equation  obtained  is  identical  in  form  viith  equation  (1) 

2 

(6)  ^  +  JH  -  m(m+l)k»^  sn^(X,k-')|y  "  0 

dx"^ 

2  2 

irtiere     k'     ■  1  -  k       is  the  complementary  modulus. 

The  characteristic  values  of  h  (and  also  H)  for  coexistence,  are 

obtained  by  equating  to  zeiro  the  appropriate  infinite  complete  denominator 

of  the  customary  continued  fraction.  When  the  continued  fraction  terminates. 

Lame  polynomials  occiir,  and  these  never  coexist,  Erdelyi  has  shown  that 

for  the  same  value  of  h,  no,  two  distinct  solutions,  one  with  2K  or  hK^ 

and  the  other  with  2iK«  and  UlK'  can  occur. 

20.   The  Lame  wave  equation 

2 
(1)    — nw-  +  (h  +  p  sn  z  +  q  an  z)w  »  0 
dz 

or,  as  it  is  sometimes  called,  the  Ellipsoidal  wave  equation  (see  [9], 

2,  pp.  159-163),  is  like  the  Herraite  elliptic  equation  in  that  it  has  four 

arbitrary  parameters.  A  physical  application  of  the  Lame  wave  equation  to 

the  problem  of  damped  electromagnetic  eigen-osciUations  was  suggested  by 

Sommerfeldo 

The  method  now  to  be  developed  applies  equally  to  the 


See  fUo],  pp.  223-22li.  Reference  is  made  to  the  work  of  M.Abraham,  [ij , 
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Ultra-ellipsoidal  wave  equation, 

(2)       — ^  +  (h  +  p  sn  z  +  q  sn  z  +  r  sn  z)w  "  0 

and  therefore  we  shall  use  equation  (2)  in  what  follows*  Let 

z  ■  am  X 
or 


/-,\  dx   /T   ,2.2  xl/2 
(3)  ^  ■  (1  -  k  sin  x)  '  . 


Hence 

/,  \    „_    k  sin  2x    „.  ^  h  ♦  p  sin  x  +  q  sin  x  +  r  sin  x    _  ^ 
(U;   w«  -  — — » K- —  w»  +  —  ■  f p    K w  ■  0  . 

aCl-k'^sin^x)  1  -  k*^  sin  x 

Now,  with  an  arbitrary  constant,  let 

(5)  w  -  u  exp(m  f     ^^?  ^^   dx) 

J   (1  -  k^sin^x) 

which  transforms  (U)  into 

U      2n 

2  m  \^^"  ^ 

/^\     _  ^  (2m-k  )3in  2x  ,  ,  o  ^ 

(6)  u»  ♦  .i 5-i — M u»  +  » M — 5"  u  =  0  • 

2(l-k^sin  x)      (1-k  sin^x) 

The  method  consists  of  dividing  (1  -  k  sin  x)   into  Y2.  ^  ^^^    ^ 

2  °   ° 

after  setting  2m  ■  k  .  The  equation  obtained  is  the  Whittaker 

equation 

(7)  u«  +  (a  +  b  sin  x  +  c  sin  x)u  ■  0 

where 

a  -  h  +  1/2 

b  -  (p-1)  -  (h+l/2)k^  -  (9kV6U)  +  2k^(h+l/2) 

c  -  r/k^ 

2 
and  also   b  -  3r  -  qk  , 

This  method  is  obviously  limited  since  the  five  parameters  of  the 
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equation  cannot  be  handled.  In  fact,  once  three  of  them  have  been  selected 
the  other  two  are  determined.  However,  it  does  provide  a  way  to  locate 
some  specific  instances  in  which  the  coexistence  question  can  be  discussed. 

For  the  Lame  wave  equation,  r  ■  0,  therefore  c  •  0,  and  hence 
equation  (7)  is  the  Mathieu  equation.  In  this  case  the  method  only 
indicates  cases  where  coexistence  can  never  occur. 
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V.    Finite  Hill  Equatjons 

21 .   Impossibility  of  coexistence  of  Ince-type  solutions 

The  solutions  of  the  Ince  equation  in  the  case  of  coexistence  are 
such  that  the  coefficients  of  the  sine  series  and  the  coefficients  of  the 
cosine  series  are  equal  except  for  at  most  the  first  N  coefficients> 
which  sometimes  all  vanish  identically.  However,  for  the  Whittaker  equation, 
as  a  result  of  the  transformation  required  to  put  it  into  the  form  of 
an  Ince  equation,  these  Ince-type  solutions  do  not  occur.  Instead,  the 
solutions  of  the  Whittaker  equation  are  Ince-type  multiplied  by  an 
exponential  in  cos  2x;.  That  Ince-type  solutions  cannot  coexist  for  a 
Finite  Hill  equation  is  the  content  of  the  follovjing  theorem. 

Theorem  21.1 

In  the  differential  equation 

(1)   y"  +  (C  +  Cpcos  2x  +  ...  *■   ^2^003  2mx)y  -  0 

where  m  >  0  is  any  integer,  two  solutions  with  period  n  of  the  form 


N  00 

(2a)  y,  -  r~  A-  cos  2nx  +  y~  A-  cos  2nx 
^   o   ^'^        H+1  2" 

and 


{2b)  y^  -  T"  B-^sin  2nx  *  Y~    A^^sin  2nx 
2   r  2n        ^1  2n 


N  00 

E 

n+l 

or  two  solutions  with  period  2n  of  the  fonn 

N  00 

(3a)     y^  ■  n  *2n-vl  °03(2n+l)x  ♦  J~    Ag^^i  co3(2n+l)x 
o  n+l 

and 

N  00 

(3^)     y2  ■  E  ^2n+l  sin(2n+l)x  +  F    A^^^^^^  sin(2n+l)x 
o  n+T 

cannot  exist  simultaneously. 
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Proof:  Asauina  that  equations   (2)   coexist  for     N  -  0.     Then  by  direct 

substitution  in  the  differential  equation  we  obtain  the  following  pair  of 

identities: 

00  P 

(Ua)     ^  £(2C  -8n  )A2j^co3  2nx  +  C2A2^cos(2n+2)x  +  ...  ♦  C2j^A2^co3{2n+2m)x 
o 

♦  C-A^  co3(2n-2)x  +  ...   +  C„  A^^cos(2n-2m)xI    -  0 

00  P 

(hb)     72  E(2C  -8n  )A     sin  2nx  ♦  C2A2^3in(2n+2)x  +  ...   +  C2j^A2j^3in(2n+2m)x 

+  C2A2^3in(2n-2)x  ♦  ..#   ♦  C2j^A2j^3in(2n-2m)x]   -  0  . 

The  recurrence  relations  obtained  from  both  (Ua')  and  (Ub)  will  be 
identical  for  n  >  m.  Hence 

(5)  ^2m*2n.2m*  —  *  ^2'^2n*2*  ^^C^"^'^^ >*2n*  ^2^2n-2*  —  *   ^2m^2n-2:n  '  ^i      ^^  >  ""^ 
Now  consider  the  recurrence  relations  for  n  ■  m. 

^^^  C2/bm*  •••  *  '^2*2m.2*  ^^""o'^^^^W   =2^2m-2*  -^  ^2m.2^2  *   "^am^o  "  ^ 

(6^^  ^2h.V  -  *  ^2^2m.2-^  ^^^'^'^^^^Zm*   ^2*2^-2*  ^"   *  ^2m-2*2  "  °  ' 

It  is  clear  that  A  ■  0. 

o 

CoriRaring  the  relationships  for  n  ■  m-1  yields  A^  ■  0,  since 
3in(-2x)  "  -  sin  2x  and  C03(-2x)  -  cos  2x.  Thus 

f"^^  C2^V-2"-*=2^2m*^2C^-3(m-l)2]A2„.2^2A2^_j^*...»C2^.^A2  -  C^^A^  -  0  . 
But  A  -  0,  and  therefore,  subtracting  (7b)  from  (7a)  gives  2C^Jl^   -  0. 
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Continuing  this  process  we  get,  upon  reaching  the  recurrence 
relationship  for  n  ■  0, 

(8)  0  -  A^  -  A^  -  Aj^  -  ...  -  A^^  . 

Then  with  (8)  and  the  recurrence  relationship  for  n  ■  1,  A^  .  ■  0 

and  30  on,  until  the  recurrence  relationship  for  n  ■  m,  A.   -  0  is  reached. 

Therefore 

(9)  0  -  A^  -  A2  -   .„  -  Ag^  -  ...  A^^  . 
Now  consider  equation  ($)  for     n  ■  ra+1 

(1°^       C2mV*2  *  ^2n.-2V  '  "'  *  °2m-2\  *  ^2«*2  '  ° 
and  applying  equation   (9) 

2m  m-*-2 
Continxilng  with     n  ■  in+2,  m+3,...  ,  it  is  clear  that  all  the     A's     must 
vanish.     Hence  the  theorem  is  true  for     N  ■  0. 

If     N  /  0,   compare  the  incurrence  relationships  for     n  ■  N^m. 
Since  by  hypothesis     A-„  /  Bj,„     but  all     A^     ■  B_       for     n  >  N*l,  we  obtain 

(11^)     C2m^2N*Um  *  —  *  ^2m-2*2N*2   "^  ^2mA2N  "  ° 

(^^^>     Sm*2N*l^  *  —   -^  ^2ni-2*2N+2  *  '^2m^2N  "  °  * 
Subtracting   (lib)  fron  (11a) 

(12)       C^(^2r  °2N>  ■  °  °^  ha  ■  ^2N 

contrary  to  hypothesis,   and  the  theorem  holds  for  solutions  with  period     k. 

If  solutions   (3)  coexist,  then  for    N  ■  0     and     n  >  «,  the  recurrence 
relations  are  the  same  for  (3a)  and   (3b),   or 

(13^     ^2m^2n.2m*l*  -  *  ^[v   i^nn)^^.^^^  ...  -  ^2^hn-2ra*l  '  ^  • 
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Compare  recurrence  relationships  for  n  ■  m-1  and  find,  by  the 
sane  reasoning  as  before,  that 

By  contlntdng  with  the  relationships  for  n  ■  in-2,  m-3,».»  we  obtain 
(lU)   0  -  A^  -  ...  -  k^^^   . 

Now,  from  the  recurrence  relations  for  n  •»  1,  we  get  A„  -  "0. 
Examination  of  the  recurrence  relations  as  n  increases  makes  it  clear 
that  all  the  A's  must  vanish. 

For  N  ^  0,  by  comparing  the  recur j-ance  relationships  for  n  ■  M  +  ra, 
we  find  the  contradiction 

(1^)   C2m(^2N.l  -  ^2Nn>  "  °  * 
which  completes  the  proof  of  the  theorem. 

22.    Impossibility  of  trigones etric  polynomial  solutions 

The  impossibility  of  the  coexistence  of  finite  trigonOTietrlc  poljrnomial 
solutions,  which  is  not  included  in  Theorem  21,''.iOuld  be  demonstrated 
by  using  the  same  procedure  as  in  Section  21,  Ho'^ever,  it  is  possible 
to  prove  the  stronger  result; 
Theorem  22.1 

Finite  Hill  equations  cannot  have  finite  trigonanetric  polynomial 
solutions . 
Proof:  Assiane  a  solution  of  equation  (21.1)  exists  in  the  form 

N  N 

(1)  y  ■  J^  ^2n°°^  ^'^       °'     ^  "  I-  ^2n'^'^  ^"^ 

where     A^„    is  the  last  non-zero  coefficient  in  either  of  the  series  (1). 
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For  n  >  m  the  recurrence  relations  ($)  holj.  Consider 
the  recurrence  relation  for  n  •  N-m,  since  all  A^^"  0  for  n  >  N+1 

(2)  Cg^A^j^  '0  or  A^j,  -  0  . 

This  contradicts  the  hypothesis  that     Apj.    is  the  last  non-zero  coefficient 
and  proAres  the  theoran  for  solutions  with  period    n. 

Assume  that  a  solution  of  equation  (2.1)  exists  in  the  form 

N  N 

(3)  y  •  E  A   3^cos(2N+l)x     or   y  »  n  A2j^^5^sin(2N+l)x  , 

o  o 

where  A^.,  ,  is  the  last  non-zero  coefficient  of  either  of  the  series  (3)« 
2N+i 

Consider  the  recurrence  relation  (13)  for  n  -  N-«n.  Since  all  A^^^^^*  0  for 
n  >  N+1 

contradicting  the  hypothesis  that  Ap  -  is  the  last  non-zero  coefficient. 
This  proves  the  theorem. 

23.    Coexistence  of  periodic  solutions 

Tn  the  discussion  of  the  Whittaker  equation  it  was  possible  to 
transform  the  Whittaker  equation  using 

y  ■  u  exp(q  cos  2x) 
into  an  equation  for  which  the  coexistence  question  could  be  completely 
settled.  The  same  basic  idea  of  transformation  enables  us  to  prove  the 
following  theoremt 
Theoran  23 ♦! 

If  in  the  Finite  Hill  equation  (21.1)  y"  ♦  (C^+CgCos  2x  ♦...♦C^jgCos  2mx)y  -  0, 
m  is  an  even  integer,  then  it  is  always  possible  to  choose  the  constants 
C  .  ...  C-_  in  such  a  manner  that  fundamental  periodic  solutions  coexist. 

O        2111 
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We  first  give  an  example  to  illustrate  the  method  more  clearly. 
Consider  the  equation 

(1)  y"     ♦   (C     ♦  C-Cos  2x  +  C.  cos  Ux  +  C/Cos  6x  +  CqCos  8x)y  "  0  . 

Let 

(2)  y  "  u  exp  j  (a  sin  2x  +  b  sin  Ux)dx 

and  obtain  the  equation 

(3)  «■  +2(3  cos  2x  +  b  cos  Ux)u'  +  (c  +  d  cos  2x  ♦  e  cos  Ux)u  ■  0 
where 

20^  -  b^  ,     C^  -  -  ab 

c  -  C^  +  a^/2  ♦  b^/2 


(U) 


d  -  C„  +  2a  +  ab 

2 
e  -  C,  -  a  /2  0 


Assume  solutions 


oo 

(5a)       ^  "  iZ.  ^?  ^^^  ^"^ 
o 


and 


00 


(5b)      u  •  T~  B2  sin  2nx 

and  get  the  following  recurrence  relations 

(e-Ub)Aj^+(d-2a)A2+2cA^  -  0 

(e-6b)A.  +  (d-ha)A, +2(c-U)A_+2dA     '0 
o  u  do 

(6) 

(e-8b)Ag*(d-6a)A^+2(c-16)A,  ♦(d+2a)A2+2eA^  -  0 

(e-10b)A^Q+(d-3a)Ag+2(c-36)A^+(d+Ua)A, +(e+2b)A2  •  0 
and  the  general  tenn  for     n  >  3 
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(7)  [e-(2n+U)b]A2^^j^+   [d-(2n+2)a>2n+2  *  ^i^'^^^^}^ 

*  [d*(2n-2)gA2^.2  *  [e  *  (2n.U)b]A2„,j.  -  0 

and  also 

(e-6b)B^+(d-Ua)B^+[2(c-l?)-(e-2b)]B2  -  0 

(8)  (e-8b)Bg+(d-6a)B^+2(c-l6)Bj^+(d+2a)B2  -  0 
(e-lJ0b)Bj^j,*(d-8a)Bp+2(c-36)B^*(d+Ua)B^*(e+2b)B2  -  0 

and   the  general  term  is  for  (n  >  3) 

(9)  [e  -  (2n*h)b]A2^^j^4[d  -  (2n*2)a]A2^^2"  2[c  -  {2n)^']^^^ 

*[A  *   (2n-2)a]A2^_2*C«  *   ^^n-UbjA^^^j.  -  0 

Now  if 

(10)  e  -  (2N+h)b    N  -  0,1,2,.., 

wo  can  set 


(U) 


*2n  *  ^2n  ^  ^    ^°^  ^  -  ^*^ 


*2n  *  ^2n  "  °    ^°^  ^  "^  ^^*^ 


and  the  following  Ince-type  solutions  with  period  n  coexist » 

00 

(12)   >  A^  cos  2nx 
H72  ^" 

00 

y  A^^sin  2nx  • 

Other  coexisting  solutions  can  be  found  for  particular  values  of  the 
parameters,  but  since  our  purpose  is  to  demonstrate  coexistence,  which 
has  been  done,  these  additional  solutions  will  not  be  elaborated. 

In  exactly  the  same  way,  assiome  the  solutions 
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00 

(13a)    u  ■  2«  Agn+l^^^^^"^*^^^ 
o 


(13b)    u  ■  H  B2j^^^sin(2n+l)x 


00 

E 

o 

and  obtain  the  recurrence  relations 


(e-5:b)A^+(d-3a)A^+[2(c-l)+(d-«)]A^  -  0 

(e-^)B^+(d-3a)B^*[2(c-l)-(d-a)]Bj^  -  0 

ilk) 

(e-7b)A^+(d-5a)A^+2(c-9)A^+[(d+a)  +  (e-3b")]A^  -  0 

(e-7b)B^+(d-5a)B^+2(c-9)B3+[(d+a)-(e-3b)jBj  •  0 

and  the  general  term  n  >  3 

(15)  t|-(2n*5)gA2^^5  *  [d.(2n*3)a]A2„^3  *   2|:c.(2n.l)2jA2^^^ 

♦  [cl+(2n-l)ajA2^.^  *  le*{2n'3)h]k^^_^   -  0 

and  the  equation  for  B  is  identical* 
Now  if 

(16)  e  -  (2N+5)b  ,     N  -  0,1,2,... 

ve  can  set 

Ag^  "^211^°       ^^^    ^  -  ^^*^ 
(17) 

k^^   ■  ^2n  "  °       ^°^    ^  "^  ^^*^ 

and  the  following  Ince-type  solutions  with  period    2n     coexist: 

00 

5~  A,  cos(2n+l)x 
n+2  "^^ 

(18)        00 

r"  A,„sin(2n+l)x  . 

This  completes  the  discussion  of  the  example. 

The  theorem  is  proved  in  exactly  the  same  way.  Make  the  transformation 
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for  the  general  Finite  Hill  equation 
(19)   y  -  u  exj)  I  (D^ 


sin  2x  +••••••  D  sin  2rx)dx 

r 


where     2r  ■  m,  and  obtain  the  equation 
(20)     Ti»  ♦  2(DTSin  2x  +  •••  ♦  D  sin  2r3c)u«  ♦  (D  +D_tCos  2x  ♦  •••  +D_  cos  2rx)u  -  0 

X  r  O      T*i.  cr 

The  general  recurrence  relationships,  with  the  usual  restriction 
n  ■  ra+r+l,  are  for  solutions  with  period    n. 


k 


(21)       [D2^.(2n+2r)Dj  "jg^^^^r         ^  ^^     ^  2[D^-(2n)2j 


L^^ 


2n 

^2n 


4  [D2j.+  (2n+2r)Dj 


3 


2n-2r 
2n-2r 


and  also  for  solutions  with  period  2n 


2n+2r+l 


(22)   [D2^.(2n+2r+l)Dj  jg*^"^^^"-^  +  2[D^-(2n*l)2j  ] 


2n+l 
^2n+l 


+  [d.  *(2n-2r>l)Dj 


■  2n-2r+l  .  q 
^2n-2r+l 


where  the  bracket  notation  indicates  that  A  or  B  can  be  used  in  the 
same  recursion  formulas.  The  criteirLa  for  coexistence  coirresponding 
to  (10)  and  (16)  are 

(23)    D2J.  -  (2n+2r)D^ 

and 
(2U)    Dg^  -  (2n*2r+l)D^  , 

and  the  corresponding  pairs  of  Ince-type  solutions  are 

00 


(25)         ^  ■  2Z    ^2n°°^  ^"^ 


n+r 


00 


u  ■  J~    Ap  sin  2nx 


n+r 
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and 

00 

(26)    u  -  ^  A2j^_^^cos(2n+l)x 
n+r 

00 

^  -  H  A  ^^^3in(2n+l)x  .  i 

n+r 

The  restriction  in  Theorem  23.1  that  m  be  an  even  integer  arises 
because  for  m  odd,  a  transformation  corresponding  to  equation  (19)  is  not 
available.  The  use  of  an  exponential  transformation  for  equations  in  which 
m  is  odd  cannot  remove  the  cos  2mx  term  from  the  coefficient  of  y 
(or  u)  without  adding  a  cos(2n+2)x  term  in  its  stead.  The  conjecture 
that,  since  m  is  odd  for  Methieu's  equation,  equations  with  ra  odd 
do  not  have  coexisting  fundamental  periodic  solutions  is  interesting,  but 
so  far  as  is  known  has  no  factual  basis. 

VI.     Singular  Equations 

2U.     Introductory  remarks 

Physical  problems  give  rise  in  a  natural  way  to  singular  equations. 
The  Schroedinger  equation,  which  describes  an  electron  in  the  electric 
field  of  a  nucleus,  for  example,  leads  directly  (see  [39],  p.2)  to  the 
singular  equation  with  periodic  coefficients 

(1)  (1-cos  2x)y"  +  (sin  2x)y'  +  (c+d  cos  2x)y  -  0. 

Equation  (l)  is  a  particular  case  of  the  Ince  equation  (11.1 )  with  a  -  b  "  1, 
and  leads  to  three-term  recurrence  relationships. 

The  investigation  of  the  oscillation  of  the  surface  of  a  lake 
(or  landlocked  sea)  by  Chrystal  led  to  his  singular  equations  which 
Halm^  -^  transformed  into 

(2)  (cos  2x)y«  +  (b  sin  2x)y«  -  (d  cos  2x)y  -  0  , 
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This  is  the  simplest  equation  with  periodic  coefficients  in  the  sense 
that  it  yields  two-term  recurrence  relationships.     The  observation  that 
both  equation  (1)  and  equation  (2)  contain  only  two  parameters  leads 
irauiediately  to  a  consideration  of  the  equation  with  three  parameters, 
which  also  has  three-term  recurrence  relationships,   namely: 

(3)       (sin  2x)y"  +   (a  ♦  b  cos  2x)y'  +  (c  sin  2x)y  -  0 

which  transfoims  into  the  hypergeometric  equation* 

The  Ince  equation  with  a  >  1  is  a  singular  equation  with 
three-term  recurrence  relationships  and  four  parameters.  By  choosing 
a  •  1,  another  equation  with  three-term  recxirrence  relationships  and 
three  parameters  is  easily  obtained: 

(U)   (1  +  cos  2x)y"  +  (b  sin  2x)y«  ♦  (c  +  d  cos  2x)y  -  0  • 

A  similar  equation  can  be  obtained  using  a  ■  -  1, 

Rellich  (see  D^lfP*^)  remarks  that  it  is  sometimes  possible  to 
change  a  singular  equation  into  a  regular  one  by  a  transformation  of  the 
dependent  or  independent  variable,  or  both;  but  he  further  states  that 
this  can  occur  in  only  a  few  simple  cases.  The  approach  adopted  here  is 
to  employ  trigonometric  transformations  to  exhibit  the  connection  between 
the  singular  equations  with  periodic  coefficients  and  well-known  equations 
in  algebraic  form.  Thus  it  will  be  shown  that  (1)  transforms  into  the 
Associated  Legendre  Equation,  (2), into  the  Dltra-spherical  or  Gegenbauer 
equation,  and  (3),  into  the  general  hypergeometric  equation.  Since 
there  exist  standard  refer ences  L-^  »  L"J  which  completely  treat  the 
solutions  to  hypergeometric  equations,  only  a  brief  discussion  of  the 
nature  and  existence  of  solutions  will  be  given* 

2$,         Two-term  recurrence  relationships.  Ultraspherical  equation* 

As  stated  above,  the  equation  (2U.2),  given  by  Halm,  leads  to 
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two-terra  recurrence  relationships,   and  consequently  may  be  considered  the 
simplest  equation  with  periodic  coefficients.     The  transformation  to  the 
Ultras pherical  equation  is  accomplished  by  letting 

(1)  t  -  sin  2x 

and  choosing     b  ■  -  Iwn     and     d  ■  Un(n+2m).     Equation   (2U.2)  then  becomes 
the  special  hypergeoraetric  equation  (see  jjo] ,  p»76) 

(2)  (l-t^)y  -   (2m+l)ty  +  n(n+2m)y  -  0 

^  IB 

whose  soltftions  are  the  ultraspherieal  polynomials   C  (t).  Therefore, 

the  solutions  of  the  Halm  equation  are  u  (sin  2x).  For  certain  values 

n 

of  m  these  become  known  orthofonal  polynomials,  in  particular,  for 
m«0,     Tchebycheff  polynoidals  of  the  first  kind 
m"l/2,   Legendre  polynomials 

m"l,     Tchebycheff  polynomials  of  the  second  kind  where  the  n  are, 
of  course,  integers. 
It  is  immediately  obvious  that  for  b"0  (i.e.,m»0),  n  need  not 

be  an  integer,  since  in  this  case  equation  (?)  reduces  to  the  trivial 

2 
periodic  equation  (setting  d"k  >  0) 

(3)  y"  ♦  k^y  -  0  . 

As  would  be  expected,  for  b  ■  -  h,   (i.e.,  m-l),  the  transformation 
(U)     y  ■  u/cos  2x 


brings  the  Halm  equation  into  the  form  of  equation  (3)  with  the  requirement 

2 

for  periodic  solutions  that  (d+U)  •  k  >  0. 

Two-term  recurrence  relationships  occur  also  in  the  equation 


See  also  Qt3j,  pp.  80-85.  Note  that  Szego  uses  a  slightly  different 
notation. 
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(5)  (sin  2x)7««  ♦  (b  cos  2x)y»  +  (d  sin  2x)y  -  0 

which  can  be  obtained  from  the  HaLn  equation  (2U.2)  by  the  phase-shift 
transformation 

(6)  X  ■  nA  -  X  . 

26.    Three-term  relationships*  The  hypergeoroetiric  and  related  equations^. 
The  transformation  of  equation  (2h.3)  into  the  hypergeornetric 

equation  is  accomplished  by  letting 

2 

(1)  t  -  cos  X 

and  by  setting 

a  »  -  UC  +  2(A+B+1) 

(2)  b  •  2(A+B) 
c  -  -  Uab 

Then  equation  (2U.3)  becomes  the  hypergeometric  equation 

(3)  t(l-t)y  ♦  [b-(A+B+l)t]y  -  ABy  -  0  . 

The  transformation 
(U)    t  -  cos  2x 

carries  equation  (2U.3)  into  the  Jacobi  equation  (see  [U3]>  p«59) 

(5)  (l-t^)y  +  [B-A-(A+B*2)t]y  +  n(n+A+B+l)y  -  0 

with 

a  -  2(A-B 

(6)  b  -  2(A+B-»-l) 

c  "  Un(n+A+B+1)  • 

The  first  of  equations  (6)  shows  the  well-known  fact  that  the 
ultraspherical  polynomials  are  Jacobi  polynomials  with  A-B. 
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The  Associated  Legendre  equation 

(7)  (l-t^)y  -  2t^  ♦  rn(n+l) 2—  ]y  -  o 

1  -  t"^ 

is  similarly  obtained  from  equation  (2U.1),  using  the  transfomaticn 

(8)  t  ■  cos  X 

with  c  ■  n(n+l)  -  2m  ,    and  d  •  -  n(n*l). 

■«. 
Equation  (2U.1)  can  be  brought  into  normal  form 


u"    +  fl/2  -  d  ♦  lA  cot  X  ♦  2i^  csc^x  u  -  0 


or 


(9)      u"     ♦  jn(n+l)  +  I  +  i  cos       +  |  +  ^  cot^^  "  sn^csc^xju  -  0     . 

2  2 

Using  the  relation  between    cot  x    and     esc  x,   equation  (9)  may  be 

rewritten  as 
(9')       u"    +    n(n+l)  ♦  r  -  (">  -  r)csc  x  u  -  0  . 

Equation  (9')  may  be  considered  as  a  limiting  form  of  Lame's 

2 
equation  for     k     ■  1« 

An  interesting  equation,   in  appearance  r«^ ambling  the  Whittaker 

equation  except  for  the  coefficient  of  the  second  derivative,   is 

(10)       (l-cos  Ux)w"   +  (h  ♦  k  cos  2x  ♦  q  cos  Ux)w  ■  0 

which  is  a  cognate  of  the  hypfergeometric  equation.  It  is  actually  the 
normal  form  of  equation  (2U.3)  which  is  obtained  b./  using  the  transformation 


/-n  \   o        /   a  ♦  b  cos  2x  ,^  ^ 
(11)   2y  -  w  exp  (J  —  gi^  ^^    d^  ; 


See  [UU],   p.  63,  where  the  phase-shifted  equivalent  of  equation  (9)  is 
given.  .  ' 
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followed  by  appropriate  modification  of  the  constants. 
Another  equation  of  the  same  type  is 

(12)  (cos  2x)w»  ♦  (h  +  k  cos  2x  +  q  cos  bc)w  «  0 
which  can  be  transformed  by 

(13)  w  «  u  exp(a  f  tan  2x  dx) 

into 
(i;*)  (cos  2x)u"  ♦  (2a  sin  2x)u'  +  (b  ♦  c  cos  2x)u  •  0  • 

It  is  not  difficult  to  illustrate  the  remark  of  Rellich  to  the 
effect  that  sometimes  singular  equations  can  be  transformed  into 
non-singular  equations.  For  example,  the  equations 

u"  *   {2   cot  x)u»  +  (c  +  d  cos  2x)u  -0 

(15)  u«  +  (b  sin  2x  +  2  cot  x)u»  +  (c  +  d  cos  2x)u  ■  0 

(1  +  a  cos  2x)u"+  (b  sin  2x  +  2  cot  x)  +  (c  +  d  cos  2x)u  ■  0 

under  the  transformation 

(16)  y  ■  u  sin  x 

become  the  Mathleu,  the  transformed  Whittaker  and  the  Ince  equations, 
respectively.  A  generalization  of  this  is  the  transformation 

(17)  y  ■  u  sin  x  cos  x 


which  transforms  the  singular  equation 

(18)      u"    *  \iL^2!LL.  .  2LJilLi]u.  ♦ 
^     '  L    sii^  ^  cos  X  J 


.(r-i)  S2^  *  s(s-l)  ^  -  r(s.l) 
sin  X  cos  X 


-  3  (r+1)  +  h  ♦  5^(x) 


u  -  0 


into 
(19)       y"+  [h  ♦  gf(x)]y  -  0 
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27 •         Four-term  reciirrence  relationships*     Description  of  the  method. 

The  method  used  to  discuss  the  conditions  for  coexistence  for 
equations  with  three-term  recurrence  relationships  consisted  of  finding 
a  recurring  2x2  determinant  which  related  the  coefficients  of  the  two 
assumed  coexisting  solutions.  As  an  indication  of  how  this  idea  might 
be  extended,  consider  the  equation 

(l)    (sin  3x  +  sin  x)w*  +  (b  cos  3x  +  c  cos  x)w'  +  (d  sin  3x  +  e  sin  x)w  •  0  . 

With  the  condition 


(2)    3b  -  c  +  d  +  e  -  8  , 


it  is  possible  to  write  the  recuirrence  relations  for  the  solutions 


(3a)   w  -  r"  A^^cos  2nx 
o   2n 


Bgj^sin  2nx 


(3b)   w  -  T^ 

with  n  sufficiently  large,  as 
(ha)   H^(A2^.  A^^.^)  .  I^(A2^.2.  A^^^j^)  .  J^U^^^,^.  A^^.,)  -  0 

(Ub)   H„(B2^.  B2„.2)  *   I„ (32^.2-^  B^^.^)  *  J^(B..^^^^*  B^^.^)  -  0 


or 


(5)   H 


n 


^2n**2n-2     ^2n*^2n-2 


*2n-2**2n-U   ^2n-2*^2n-U 


n 


^^n-U*^2n-6   ^2n-l4*^2n-6 


^2n-2*^2n-U  '   ^'2n-2**2n-U 


Multiplying  ($)   by  (A2j^.2**2n-U^^^2n-2'*'^2n-U^*  ^"^  rewriting  in 
determinant  form,  we  get 


(6) 


A„  +. 


2n*^n-2  *2n-2*^2n-U 


^2n*^2n.2  ^2n-2*®2n-li 


J 
n 

w 

n 


^2n-2*^2n-U  *2n-U**2n-6 


^2n-2*^2n-U  ^2n-U*^2n-6 


which  can  be  written  as  3x3  determinants 
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(7) 


*2n  ^2n-2  ^2n-U 


^2n  ^2n-2  ®2n-ii 


-1 


J 
n 


-^2n-2  *2n-U  *2n-6 


^2n-2  ^2n-U  ^2n-6 


-1 


where  H  -  Un  -  2bn  -  d 


and 
(8) 


J^  -  -  U(n-3)  -  2b(n-3)  ♦  d  . 


A  similar  analysis  could  be  developed  for  solutions  with  period  n» 
If  it  were  known  that  the  equation  had  periodic  solutions  with 

period  n,  equation  (?)  could  be  used  to  examine  the  conditions  for 

coexistence. 

28 •    Remarks  on  the  periodic  solutions* 

The  statement  that  a  differential  equation  with  singularities  has 
periodic  solutions  is  meaningful  only  if  the  periodic  solution  is 
continuous  everywhere  and  has  a  continuous  first  derivative,  or  if  the 
periodic  solution  is  an  analytic  function.  Equation  (2U.3),  which  is 
equivalent  to  the  hypergeometric  equation, will  have  two  periodic 
solutions  (in  general  with  singularities)  if  the  solutions  are  single- 
valued  analytic  functions.  For  this  it  is  necessaxy  but  not  sufficient 
for  the  exponents  of  the  hsrpergeometric  equation  to  be  integers  or  half 
integers.  Even  if  the  exponents  are  integers,  the  solutions  may  be 
logarithmic,  and  periodic  solutions  will  not  occur  (see  [9],  pp.  71-73)« 

In  his  discussion  of  the  frequency-modulation  equation, 
(1  ♦  a  cos  2x)y»  ♦  cy  -  0 
Cambi'^-'  points  out  that  for  a  >  1,  this  becanes  the  Legendre  equation. 
For  a  ■  1,  he  states  that  ajl  solutions  are  stable  and  of  period  2k» 
However,  he  does  not  show  that  there  is  even  one  solution  regular  at  both 
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singTilarities,  although  it  is  a  curious  fact  that  the  algebraic  form  of 
the  equation  aoes  have  two  linearly  independent  solutions  which  are  finite 
at  both  finite  singularities.  As  Cambi  indicates,  the  solutions  have, 
in  general,  first  derivatives  with  logarithmic  singularities.  Thus  the 
solutions  are  not  periodic  in  the  usual  sense. 

VII.   Conclusion 

29 •         Summary  of  results. 

The  coexistence  problem  for  Hill's  equation,   that  is,   the  problem 
of  establishing  the  conditions  under  which  two  independent  solutions  of 
period     n     or  period     2n     exist  simultaneously  for  the  equation 

y"  +   [F>  +  gl(xT\j  -  0 
where     (?(x)     is  of  period     n,  has  been  investigated,   using  recurring 
2  by  2     determinants.     The  determinants  are  formed  bv  equating  the 
recurrence  formulas  of  the  coexisting  solutions,   and  the  criteria  for 
coexistence  are  conditions  on  the  parameters  of  the  equation,   determined 
from  these  recurring  determinants.     The  method  is  applicable  to  equations 
from  which  three-term  recurrence  relationships  can  be  obtained  either 
directly  or  by  transformation  of  the  dependent  or  independent  variable. 

The  coexistence  problem  is  intimately  connected  with  the  stability 
or  instability  of  the  solutions  of  the  equation.     The  solutions  of  period 
n     and     2n     provide  the  boundary  curves  between  the  intervals  of  stability 
and  the  intervals  of  instability  in  the  stability  diagram.     The  coexistence 
of  any  two  of  these  fundamental  periodic  solutions  implies  the  vanishing 
of  the  correspondiiig  interval  of  instability. 

The  previously-known  results  on  coexistence  for  the  f'athieu  equation, 
the  Lame  equation,   p.nd  the  Frequency-modulation  equation  are  shown  to  be 
special  cases  of  the  solution  of  the  coexistence  problem  for  the  new 
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equation  (with  |a|  <  1) 

(1  ♦  a  cos  2x)y"  +  (b  sin  2x)y'  ♦  (c  ♦  d  cos  2x)y  ■  0 
irtiich  we  have  named  the  Ince  equation,  in  honor  of  E.  Lindsay  Ince,  the 
British  mathematician.  The  Ince  equation  has  four  arbitrary  parameters, 
and  is  thus  more  general  than  the  previously  treated  equations  which  had 
only  two  or  three  arbitrary  parameters.  The  Ince  equation  is  the  most 
general  non-singular  equation  with  trigonometric  coefficients  which  leads 
directly  without  transformation  to  three-term  recurrence  relationships. 
It  is  shown  that  certain  equations  with  elliptic  functions  as  their 
periodic  coefficients,  such  as  the  Lame  equation,  and  the  Hermite  elliptic 
equation,  can  be  transfonned,  by  a  change  of  independent  variable,  into 
an  Ince  equation. 

It  is  proved  that  for  coexistence  to  occur  for  the  Ince  equation, 
it  is  necessary  and  siifficient  that 

p+q-*n   or  -  1/2         m,n  ■  0,1,2,  ... 

2  1/2 

where     p  -  bA^     and     q  ■   (b     ♦  had)  '    Aa 

and  c  has  an  appropriate  characteristic  value. 

There  are,  in  general,  three  tjrpes  of  coexisting  periodic  solutions 
which  may  occur:  truncated  solutions  for  which  the  first  N  coefficients 
of  the  odd  and  even  infinite  trigonometric  series  vanish,  split  solutions 
in  which  the  first  N  coefficients  are  different  for  each  of  the 
trigonometric  series,  and  polynomial  solutions  whi  ch  may  be  either  truncated 
or  split. 

If  p  -  q  -  n,  even  intervals  of  instability  will  vanish,  while  if 
p  -  q  B  m-1/2,  odd  intervals  of  instability  will  vanish.  Both  even  and 
odd  intervals  of  instability  can  vanish  leaving  only  a  finite  number  of 
intervals  of  instability  if  p  -  q  ■  2n  and   m  -  1/2,  simultaneously. 
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The  simplest  example  of  this  is  the  equation 

(1  +  a  cos  2x)y"    +  (b  sin  2x)y'  +  cy  -  0 
for  which  all  even  intervals  of  instability  (except  the  aero-th)  vanish, 
and  for  which  odd  intervals  can  vanish  vrtienever 
b  -  -   (2n*l)a  , 

The  Frequency-modulation  equation 
(1  ♦  a  cos  2x)y''  +  cy  -  0 
is  an  example  for  which  all  even  intervals   (except  the  zero-th)  vanish,  but 
for  which  the  odd  intervals  can  never  vanish. 

These  examples  illustrate  the  fact  that  all  even  intervals  of 
instability  (except  the  zero-th)  vanish  for    d  -  0,     Similarly  for 
d  ■  a+b,  all  odd  intervals  of  instability  vanish,  while  for     d  •  0     and 
a  ■  -  b,  all  intervals  of  instability  (except  the  aero-th)  vanish.     This 
last  case  corresponds  to  a  Hill  equation  with    (S^(x)     a  constant. 

If  the  coexisting  solutions  obtained  for  the  Ince  equation  are 

designated  Ince-type  solutions,   it  is  pi*oven  that  Ince-type  solutions 

cannot  coexist  for  the  Finite  Hill  equations 

M 
y"  *(  2^  ^2m°°^  2mx)y  -  0 
o 

and  also  that  trigonometric  polynomial  solutions  cannot  occur. 

If  m  is  an  even  integer,  it  is  shown  that  coexistence  of  fundamental 
periodic  solutions  can  occur.  The  solutions  are  Ince-type  multiplied 
by  an  exponential  function  with  a  periodic  exponent. 

The  particular  Finite  Hill  equation  with  ra  «  2  is  called  the 
Whittaker  equation,  and  is  the  next  in  increasing  complexity  in  the 
sequence  of  Finite  Hill  equations  after  the  Mathieu  equation.  The  Whittaker 
equation  is  most  conveniently  written 

y»  +  (h  ♦  bnq  cos  2x  •••  2q  cos  Ux)y  ■  0  , 
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The  condition  for  coexistence  is  that  m  be  a  positive  or  negative 
integer  or  aero,  with  h  the  required  characteristic  value.  For  q 
imaginary,  the  condition  for  coexistence  is  unchanged.  If  ra  is  an  even 
integer,  odd  intervals  of  instability  will  vanish,  and  conversely,  if  m 
is  an  odd  integer,  even  intervals  of  instability  will  vanish. 

A  number  of  singular  equations  were  examined  and  identified  with 
known  equations,  such  as  the  hypergeometric  equation.  The  idea  of  using 
2  by  2  determinants  was  extended  to  a  method  of  employing  3  by  3 
determinants,  wl-dch  can  be  used  to  treat  scsne  particular  singular 
equations  leading  to  four-term  recurrence  relationships. 
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The  Glaeslflcatlon  of  Linear  Differential  Equations 
Linear  second  order  differential  equations  with  rational  coefficients 
can  be  classified  in  a  systematic  way  according  to  the  kind  and  number 
of  singular  points  of  the  equations  (see  [26],  p.  U9U).  This  scheme  stems 
from  the  Klein-Bocher  discovery L  -J  in  18 9U  that  many  important  problems 
in  mathematical  physics  lead  to  equations  which  can  be  derived  from 
a  differential  equation  with  five  elementary  singularities.   Six 
equations  can  be  derived  from  the  equations  with  five  elementary  singulari- 
ties by  a  suitable  coalescence  of  two  or  more  singularities.  These 
derived  equations  are  the  equations  of  Lame,  Mathieu,  Legendre,  Bessel, 
Weber(or  Herraite)  and  Stokes  (a  particular  case  of  Bessel's  equation). 
Ten  equations  can  be  derived  from  the  equation  with  six  elementary 
singularities,  fourteen  equations  from  the  equation  with  seven  singularities; 
and  this  may  be  continued  for  any  number  of  elementary  singularities. 
There  is  no  loss  of  generality  if  three  of  the  singularities  are  fixed  at 
the  points,  0,  1  and  oo,  since  any  three  points  can  be  transformed, 
without  change  in  the  exponents  relative  to  these  points,  by  a  projective 
transformation  of  the  independent  variable. 

The  method  of  classification  developed  here  is  based  upon  the 
Klein-Bocher  discovery  and  upon  the  observation  that  differential  equations 


An  elementary  singularity  is  a  regular  singularity  with  exponents 
(relative  to  the  singular  point)  which  differ  by  1/2,  The 
coalescence  of  two  elementary  singularities  gives  rise  to  a 
regular  singularity  with  an  arbitrary  exponent  difference,  while  the 
coalescence  of  three  or  more  yields  an  irregular  singularity. 
Irregular  singularities  are  further  characterized  by  their  specie| 
an  irregular  singularity  of  the  m-th  specie  is  one  generated  by 
m+2  elementary  singularities. 
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with  periodic  coefficients  can  be  transformed  into  algebraic  differential 

2 
equations  by  a  trigonometric  trasnformation  (e»g.,  tscos  x)  or  by  an 

2 
elliptic  transforoation  (e.g.,  t«sn  x).  These  algebraic  equations  can 

be  derived  (by  coalescence  of  singularities)  from  differential  equations 
with  elementary  singularities.  The  purpose  of  this  new  classification 
method  is  to  display  the  differential  equations  with  periodic  coefficients 
in  a  manner  which  shows  their  relation  to  one  another  and  at  the  same 
time  indicates  their  position  in  the  general  scheme  of  linear  second 
order  differential  equations. 

To  accomplish  this  purpose  a  sequence  of  three  dimensional 
manifolds  is  employed.  The  accompanying  charts  show  four  manifolds 
of  interest,  the  Spherical,  Spheroidal,  Ellipsoidal,  and  Hyper- 
ellipsoidal.  These  manifolds  are  best  characterized  and  more  easily 
understood  by  referring  to  the  general  differential  equations  corresponding 
to  each  manifold.  These  general  differential  equations  are  given 
in  both  algebraic  and  trigonometric  form  in  Table  I.  Since  three 
singularities  may  be  fixed  without  loss  of  generality,  three  of  the 
singularities  are  located  at  the  points  0,  1,  and  oo.  The  equations 
in  the  spherical  and  spheroidal  manifolds  have  three  singularities  , 
irtiile  those  in  the  ellipsoidal  and  hyper-ellipsoidal  manifolds  have  four 
singularities.   This  fourth  singularity  may  be  located  at  an  arbitrary 
point  and  is  a  parameter  or  irreducible  constant  of  the  equations.  The 
three  axes  used  in  the  manifolds  are  §,  r  and  nj  where  0  is  twice 
the  exponent  difference  at  the  origin;  r  is  twice  the  exponent  difference 


These  singularities  may  be  elementary,  regular  and  irregular.  It  shoulQ 
be  noted  that  the  equations  in  a  given  manifold  are  not  all  derivable 
from  the  same  equation  with  only  elementary  singularities. 
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at  1;  and  n  is  the  index  in  the  coefficient  of  w.  The  specie  of  an 
irregular  singularity  at  infinity  is  related  to  n. 

Individual  equations  are  designated  by  a  capital  letter  and  three 
subscripts,  9,  r,  and  n.  T  is  used  for  the  spherical  manifold  and  S 
for  the  spheroidal,  E  and  H  are  used  for  the  ellipsoidal  and  hyper- 
ellipsoidal  manifolds,  respectively,  when  the  fourth  singularity  is 
regular;  and  L  and  N,  respectively,  when  the  fourth  singularity 
is  elementary.  In  the  acccxnpanying  chairts,  many  of  the  equations  are 
named.  In  general  these  correspond  to  the  names  given  in  the  literature, 
although  a  few  appear  for  the  first  time.  Equations  Lq,,  and  L^  ^ 
were  given  by  Hermite  in  1880  but  are  not  yet  named.  In  the  case  of  the 
hyper-ellipsoidal  manifold  only  one  equation  appears  in  the  literature, 
%12  '"•nfcioned  by  Ince  in  1923.  The  names  indicated  for  N  and  H 
equations  should  be  considered  tentative.  T,  _  was  named  by  analogy  with 
T-  -,  the  Associated  Mathieu  equation,  but  as  may  be  noted,  is  not  in 

consonance  with  S,--,  the  Associated  Legendre  equation.  Actually  S^^, 

2 

must  be  transformed  by  t»a   into  the  Associated  Legendre  equation, 

but  the  connection  seems  close  enough. 

From  the  standpoint  of  differential  equations  with  periodic 

coefficients  the  sequences,  T,,  ,  L,^  ,  and  E..^   appear  to  be  most 

important.  It  is  interesting  to  observe  that  the  coexisterice  problem  is 

solved  completely  in  all  three  cases  for  n-1,  but  for  n"2,  is  solved  only 

in  the  case  of  T.,o»  Also,  except  for  some  general  statements  (Borg's 

theorems)  and  our  results  for  T,,  with  n  even,  little  is  known  for 

lln  * 

arbitrary    n.     In  terms  of  the  Kiein-Bocher  classification,  equations 
T^^^     (Mathieu »s  equation)  and    L-, ,      (Lame's  equation)  may  be  derived 
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from  the  equation  with  five  elementary  singularities.  Equations  T^^^ 
aivl  E,  ,  are  derivable  from  the  equation  with  six  elementary- 
singularities.  Equation  T^,^  ^^^  ^^°   elementary  singularities  and  one 
irregular  singularity  of  the  n-th  specie  and  therefore  may  be  derived  from 
the  equation  with  n+U  elementary  singularities. 

The  major  new  results  of  this  paper  concerned  equations  T^-jp* 

E,,,   and  T,,   (n  even).  The  classification  scheme  developed  here 
111        lln 

would  have  served  a  useful  purpose  merely  by  exhibiting  the  relation  of 
T, ^P  and  T^.   to  the  hypergeometric  equation,  and  of  B^,  to  Heun's 
equation.  By  displaying  all  of  the  equations  discussed  in  a  systematic 
way,  the  classification  seems  to  provide  a  coherence  to  the  discussions. 
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CUSSIFICATION  OF  EQUATIONS     -     SPHERICAL  MANIFOLD 
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Chart  II 

CLASSIFICATION  OF  EQUATIONS  -  ELLIPSOIDAL  MANGOLD 
(Exponent  difference  of  singularity  at  a  equal  to  1/2 ) 
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Table  I  -  Classification  of  Equations  -  Algebraic  and  Transcendental  Forms 


SPHERICAL  -  Tq^^: 


^Af 


w  ■  0 


.   Pd    r  Iw   o"  " 

w"  +  [(r-l)cotx  -  (9-l)tanx]w«  -  ij^  A^cos^^x  w  -  0 


SPHEROIDAL  -  S-,  : 
— — — — —   9rn 
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ut(t-ir 
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sin  X 


ELLIPSOIDAL  -  Eq^: 
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HYPER-ELLIPSOIDAL  -  H, 
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